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RELATIONSHIP BETWEEN ION DISTRIBUTION AND MEMBRANE 
POTENTIAL DURING A STEADY STATE 


DANIEL L. GILBERT* 
DEPARTMENT OF PHYSIOLOGY 
ALBANY MEDICAL COLLEGE OF UNION UNIVERSITY 
ALBANY, NEW YORK 


Equations were derived showing the relationship between the membrane 
potential and the quantities which influence it under steady state con- 
ditions. Essentially, the membrane potential is caused by the valence 
and concentration of the non-permeating ions. The permeating ions can 
modify the membrane potential by altering the relative concentration of 
the non-permeating ions with respect to the concentration of the permeat- 
ing ions. 

For muscle, the sodium cations act as the non-permeating ions in the 
extracellular environment by the maintenance of some type of active 
metabolic process and large anions act as the non-permeating ions in the 
intracellular environment. Both of these non-permeating ions contribute 
about equally to the maintenance of the resting membrane potential. When 
the active metabolic process for sodium extrusion breaks down or when 
acids are added, the membrane potential should decrease. Water should 
enter the cell when the sodium metabolic process is diminished; water 
should leave the cell when acids are added. When acid is added, it is 
expected that the cations potassium and sodium will leave the cell with 
little or no shift of the chloride ions. 


Introduction 


It has been observed that the ionic constituents within the bio- 
logical cell are different from the ionic constituents of the cellular 
environment. It has also been observed that there is an electrical 
potential existing across the cell membrane. One of the factors 
responsible for both this ionic difference and the membrane po- 
tential is the presence within the system of non-permeating ions, 
It is the purpose of this paper to describe mathematically a system 


*Present address: Department of Physiology, Jefferson Medical Col- 
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which contains both permeating and non-permeating ions during a 
steady state, and also to show the relationship of these ions to the 
membrane potential. In the first part of this paper, the general equa- 
tions describing such a system are derived. The second part de- 
scribes some of the effects of changing the valence of the non- 
permeating ions upon the system and also illustrates how these re- 
lationships can be applied to muscle physiology. 


General Equations 


Dependence of a given permeating tonic species upon the mem- 
brane potential. The relationship between the flux of a given species 
and its potential gradient is linear (Onsager and Fuoss, 1932) and 
can be expressed as: 


Ui eases (1) 


a) y = flux of a given species 

b)C = concentration of the species (for rigorous treatment, 
activity should be used instead of concentration) 

c) B = diffusion coefficient of the species 

d) PR = gas constant 

e) T = absolute temperature 

f) s = distance 

g) w = total potential of species 


The total potential @ can be defined as: 


w=u,t+RTInC+2FE,+W (2) 
a) wu, = Standard chemical potential of species 
b) 2= valence of species 
c) F =the Faraday 
d) E = electrical potential 


e) W = potential of species due to other forces besides elec- 
trical and chemical ones. 


When some barrier (membrane) separates the internal environ- 
ment (compartment 1 or side 1) and the external environment (com- 
partment 2 or side 2), then: 


viet ws (3) 
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It 


a) y, = the net flux of species going across the barrier from 
compartment 2 to compartment 1 
b) y, = the flux of the species going from side 2 to side 1 


C) y, = the flux of the species going from side 1 to side 2 


Levi and Ussing (1948) have derived a relationship between flux 


; ae dE 
and concentration. In their derivations, they assumed that —¥ was 


constant and did not include the W term. If the W term is included 
dW 

and if it is assumed that on is also constant, then it is possible to 
8 


determine the more general relationship between flux and concen- 
tration, uSing an analogous derivation. For mathematical conven- 
ience, some other symbols are used which are defined by the fol- 
lowing equations: 


E=E,,-£,; (4) 


a) Subscript 1 refers to compartment 1 
b) Subscript 2 refers to compartment 2 
c) FE = membrane potential 


Wo =W,-W, (5) 


a) Wo= work potential of species. In biological systems, W, 
can refer to the work derived from metabolism. 


RT 
E = —Inr (6) 
salt 
W,=RT lng (7) 
K = of (8) 
on 
a) ¢ = thickness of membrane 
b) K = permeability coefficient of individual species 
og inne 
spy tome (9) 
r?qg—1 
Incp* 
k= saestalid (10) 


r7qg—1 
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The general relationship between flux and concentration is: 


Y,= yl, -h,C, =k, (r79C, - Cy) (11) 
aly, =k,C, 
b)y, =%,C, 
The relationship between current density and flux is: 
J, = aFy; (12) 


a) J, = net current density of individual species going from 
compartment 2 to compartment 1. 


For the existence of a steady electrical state, the summation of all 
the net current densities must be equal to zero. For the existence 
of a steady chemical state of a given ionic species (in which there 
is no production or consumption of the ionic species), the net flux 
of this species is equal to zero. If there exists a steady chemical 
state for all the ionic species, then a steady electrical state must 
also exist. However, a steady electrical state can exist in a sys- 
tem without the existence of a steady chemical state. The Hodgkin- 
Katz (1949) equation showing the relationship between permeating 
ionic species and the membrane potential is based upon a steady 
electrical state, and therefore for steady chemical states, their 
equation is not useful (Conway, 1957). If equation (11) is set 
equal to zero, then a steady chemical state exists which results in 
the following (Gilbert and Fenn, 1957): 


oe = r*q (13) 


Equation (18) shows the dependence of a given ionic species upon 
the membrane potential (r is a function of the membrane potential, 
see equation (6)), but does not interpret the cause of the membrane 
potential. If W, is equal to zero; then from equation (7) it follows 
that g = 1; and then (13) describes the familiar Donnan equilibrium 
(Donnan, 1924). If the activity coefficient of the ionic species is 
the same in both compartments, then the use of concentration units 
in (13) will not introduce error. However, if the activity coeffi- 
cient is not the same, instead of using only g in the equation one 
should use q times the ratio of the activity coefficient for compart- 
ment 2 to the activity coefficient for compartment 1. 

Relationship between redox potentials and membrane potentials. 
Since all solutions possess redox potentials, it is desirable to 
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Show their relationship to the membrane potential. The electrical 
potential observed between the two compartments is given by the 
following: 


B=Ek +E. (14) 
a) E = observed electrical potential 
b) BE, =E, - E,,, = reduction-oxidation potential 


If a solution is at equilibrium, then the reduction-oxidation potential 
of each species will be equal to the observed reduction-oxidation 
potential of the solution. The value of the reduction-oxidation po- 
tential of a given species is: 


(15) 


a) &,,, = reduction-oxidation potential of individual species 

b) BS ay = standard reduction-oxidation potential of individual 
species 

c) oxid = oxidized form of species 

d) red =reduced form of species 

e) m = number of electrons liberated in the transformation 
of red to oxid 


Sometimes, biological workers use another nomenclature: 


EB, =#,,-#,. =- Ey = £2. -=,,,1 = oxidation-reduction po- 
tential (16) 


At equilibrium, EF = 0, so from equation (14): 
a Soe Oe (17) 


This equation (17) has been pointed out to exist for permeating 
ions which can undergo oxidation-reduction (Korvezee, Dingemans, 
and Dykgraaf, 1947). 

General conditions necessary for a steady-state. Again, it should 
be pointed out that the equations (6), (13), and (17) do not give any 
interpretation as to the cause of the membrane potential, but merely 
point out relationships which must exist. One of the main pur- 
poses of this paper is to derive an equation which will give a 
meaning to the membrane potential during a steady state. 

Various conditions of restraint must be considered for the steady 
state. When there is no net movement of water; the following must 


be true: 
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D,+2,=D,+2,=S8 (18) 


a) D = concentration of permeating substances 
b) aw = concentration of effective non-permeating substances 
c) § = total concentration 


The effective non-permeating substance concentration is composed 
of: 
e=X+L (19) 


a) X = actual concentration of non-permeating substances 
The value of L is defined as: 
L =— 20 
RT C9) 


a) ]=net pressure on opposite side forcing water into same 
side, 


For osmotic equilibrium, /= L =0, so that =X. For simplicity, 
only osmotic equilibrium conditions will be considered. Let the 
value of g be defined as: 


“1 (21) 


ae 
X, 


Equation (22) takes into account the water volumes: 
Vie te a Vee ee (22) 


a) V = volume of water 
b) subscript @ refers to new steady state after some change 
c) subscript e refers to increase in quantity 


By definition: 
f= V 23 
a ( 


If f is less than one, water was added to that side; if f is greater 


than one, water was removed from that side. If fg =lande,, =0 
then: * 


la Ga 
g 


(24) 
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a) X., = amount of non-permeating substance added to side 1 
times the initial volume of side 1 

b) X,, = amount of non-permeating substance added to side 2 
times the initial volume of side 2 


If the diffusible substance is a uni-univalent electrolyte and if 
g = 1, then equation (13) leads to the Donnan equilibrium: 
M A 


+= 2 (25) 


r= — 
M, 1 


a) M refers to the concentration of permeating univalent cations 
b) A refers to the concentration of permeating univalent anions 


Another ratio can be obtained from equation (25): 


M M 
fae (26) 
A, A, 
Electrical neutrality must also be taken into account: 
M-A+nX =0 (27) 


a) n = valence of X. Since X is composed of all the non-per- 
meating substances, n is an integrated mean valence. 


The concentration of the non-permeating substances in equivalents 
per unit volume is given by @ or P as follows: 


G=-P=nX (28) 
The summary of the general conditions of equilibrium are: 
M,+A,+%X,=M,+4,+%,=8 (29) 
M,A, =M,A, (30) 
A, =M,+7,%, (31) 
A, =M, +7 %X, (32) 


From the summary of the general conditions for equilibrium, the 
following relationships can be derived: 


a M, My My Ign, +1) -41? - 4," (33) 
oa MeocA, A, = [g(n, - 1) + 112 - 0,2 
2 2 
cg ects Seas (34) 


AS 4(1- 9) 
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1 2 
—(n +y-1 See 
oa Eps bee (35) 
: Re -1)+ | -1, 
g 
2 
Fo. +1] -n, 
aes (36) 


A. Conditions of restraint. By not permitting any concentration 
to have a negative value, it is possible to show from equations 
(33), (34), (35), and (36) that one of the four following conditions 
must occur: 


Condition lis: -1<g-1<0<|n,|<|a,|+1)-1<|n,| 


1 
Condition II dave i" Oat fe = tie 41) —1 ape 
g g 


1 1 
Condition T,tas-1.<-—1. = O's ny b= (1p lated ie hentia 
g 


If a constituent concentration equals zero, then Condition IV 
occurs and is: 


1 
Osim, le sag 1 #2) —1 


If one constituent concentration equals zero, then at least another 
constituent must equal zero. Subconditions of Condition IV are: 
Subcondition A, Condition IV is: 0< g<1and0<| Ns |<| N, | 
Subcondition B, Condition IV is: 1 < g < « and0<|n,|<|n, | 
Subcondition C, Condition IV is: g = 1 and 0 <| n, | = [n, | 
The value of g must be greater than zero and less than infinity. 


The trivial cases of X, =X, = 0 and/or D, = D, = 9 will not nec- 
essarily be defined. 


B. Relationships between the variables, The relationships be- 
tween r and /A and n,, nq, and g are: 
P (Mg F ng) 3 (1 a 9) 


"(n, +n) +(1=9) oo 
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, - ta ~ 19) + (1-9) 


38 
(ng —n,9)- (1-9) ea 
If M, = M, =0, then: 
1 
ee 
n 
r= (39) 
1 
1+ — 
No 


This condition can occur when mM is not present even though the 
membrane is permeable to M. If a membrane is not permeable to 
any cations, then the same condition can occur. Thus, for red 
blood cells, the cations are effectively non-permeating due to 
metabolic processes and therefore equations (39) and (40) can 
apply. 

Similarly, if A, = A, = 0, then: 


1 
7 Se 
nN 
2 
f= 41 
2 (41) 
se ee: 
ths 
h— © (42) 


According to the Conway-Boyle theory (Boyle and Conway, 1941) 
of muscle permeability, the potassium and chloride ions obey the 
Donnan equilibrium. If the chloride is removed from the system 
and replaced by a non-permeating ion, then equations (41) and (42) 
apply. 

If M, = A, = 9, then: 


? —» 00 (438) 
1 
1+ — 
n 
h= 2 44 
2 (44) 
5 Rema 
nN 


When the potassium concentration in the extracellular fluid of mus- 
cle is decreased towards zero, then equations 2) and (44) are ap- 
proached as a limiting condition. 
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Similarly, if M, = A, = 0, then: 
r=0 (45) 


De 
n 
hie - 46 
: (46) 
5 Care 
Mo 


Whenever g = 1 (i.e., Condition III or Subcondition C, Condition 
IV), then the following must be true: 


(A-1)(r-1)=0 (47) 
0<|a,| = In, | (48) 
de. = Me > 0 (unless r or h — o, and then 
1-rh r—-h ° 
ny No 
1-rh r—-h Mets? 


If g is not equal to 1, then r and/or / are only functions of n,, 
m,, and g. Thus, from equations (6) and (37), it can be seen that 
the membrane potential is a function only of the valence of the 
non-permeating ions and the ratio of the concentration of the non- 
permeating ions in one compartment as compared to the other (as- 
suming that the temperature is constant). If there is a non-permeat- 
ing species present on one side, there must necessarily be a non- 
permeating species on the other side. It can be seen from the 
conditions of restraint that the value of g can only be within a 
given range depending upon the valences of the non-permeating 
species. If the valence of the non-permeating species is zero 
(i.e., n, =n, = 0), then from equations (47), (48), and (49), it can 
be seen that r = h = 1; thus the membrane potential is equal to zero 
(equation (6)). The membrane potential is therefore due to the 
valence of the non-permeating ionic species and so equations (6) 
and (37) give an interpretation as to the cause of the membrane po- 
tential. When g=1, then r and/or A is not restricted by just n, 
and n, and the above equations are not useful. However, if some 
additional mathematical relationships are made, then these equa- 


tions can be used even when g is equal to one (excluding the 
Special case when n, =n, = 0). 
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M m(a-1)+1 a(6-1)-1 
n(a-1)-1 n,(b-1)+1 


(50) 


_n(e-1)-1 n(d-1)-1 
Fe(Oes amene ld 1). 4 1 


The symbols a, 6, c, and d are related to M1» N,, and g as follows: 


1 n 
(1-=) a= (1+ 2) (52) 
g Ne 


(51) 


Org) tn [4 = (53) 
ill ny 

(t- 2) e+ (1-22) (54) 

GQaga- (1-22) (55) 


The values of a, b, c, and d represent the ratio of n, to n, and g 
to one, even when g is equal to one. 

As pointed out in the above relationships, the concentration of 
the non-permeating ionic species is the cause of the membrane po- 
tential. It would therefore be expected that the greater the con- 
centration of the non-permeating ionic species (with the same 
valence), the greater should be its influence upon the membrane 
potential. The magnitude of the concentration of the non-permeating 
ionic species has more meaning if used in comparison to the con- 
centration of the permeating ions. Therefore, equations can also 
be derived showing the relationships of 7, h, and g to n, andn, 


Xx 


iy ees 1 
and the concentration ratios, Fis and 5s 


2 2 


X = 
Qn, — n(n, - 1) The -Vo 
2 


oo Sg es one 
2(n, - 1) gaa 
2 


2(n, + 1) 


r ze 
2n, — n(n, - 1) Te +Vo 
2 


f= 


(56) 
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X 
2n, — n(n, - 1) Fe -Vo 


h= : 
2(n, - 1) 


2( 1) ( | 
Mig oe =<. Sa 
: 1 2 A, 


X 
2n, -7,(n, - 1) 7 +V@ 
2 


X 7s X 
2—(n, -1)—*-YVo rary ¢ Seas) Ge 
A, A, 
2 (ie aes gee Te eee 
1 A, 2 A, 
where 
Aa X 
a=," (ng — 1)? —* + dn? —n,) 2 +4 
Ay 
xX, 
an, +n,(n, - 1) 7° + Ve 
Ppp teen a Ba lle rt tae 
Xx, 
2(n, + 1) + 2n,(n, - 1) — 
A, 
2(n, - 1) 
“& X, 
ang + ni(ny - > - ve 
2 
Xx, 
an, = Wi Ve 
| en ee: A 
2(n, — 1) 
Xx, 
2(n, + 1) - 2n,(n, + 1) — 
A, 
Xx, 
an, ~n,(ng 21) Sve 
2 
x Xx 
2+(n,-1) 7° + Ve (n, — 1)? —+ 
rer 2 
44(n?2-1 Xx, 9 Xx, 
os Te i - og Day = 
2 2 


(57) 


(58) 


(59) 


(60) 


(61) 


(62) 
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where 


2 o Ay x, 
€=n,(n, - 1) het (note) a (63) 
A, A, 

It is evident from equations (56) and (59) and also from equa- 
tions (60) and (63) that when the concentration of the non-permeating 
ionic species approaches zero (or conversely when the concentra- 
tion of the permeating ionic species approaches infinity); then the 
value of 7 approaches 1 and thus the membrane potential approaches 
zero, aS expected, 

A determination of the concentrations on one side and the mem- 
brane potential permits the calculation of the net valences of the 
non-permeating ionic species by solving for n, in equation (32) and 
n, in equation (56): 


ON eaters (64) 
2 2 
[ns 32- | p41 
= 9 (65) 
he A, 


Steady state conditions not assuming ‘‘g’’ is equal to one. If 


there is no permeating ion which obeys the Donnan equilibrium, 
then equation (25) is not valid, and instead the following equations 
and definitions can be used: 


M 
i, =1dy = ly (66) 
nS a r, (67) 
A, 
Let 
Q= Wuda (68) 
M.A, =M,4,@ (69) 
A 
ast ge ae ay (70) 
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The relationships between the variables are given by the fol- 
lowing equations: 


[e3)- Ci orsb td 
[eet) eae 


Xx 
(n, - 1) (1 n, a Qr 4” = g -» 2-9] Rit, 
2 


2 
+(n,+1)=0 (73) 


2 


+ A(1 ~ Q) (i= cal [ies - 0932 


; ie 0 (74) 
2 A, 


~ [2,0 +n,(n, - 1) Z| r,+(ng-1)=0 (75) 
2 


X iar 
saa ean 25] 9 


2 


Ane A 
+2 | 2n(@-1) 42+ (ay -1) ibe 


ge 
+ (n, - 1) ree (76) 
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2 


[Ig(in, - 1) + 11? -n,7} Q tobe 
A? 
1 


+ 2[971Q(n,? —i)- 2n 7] + [29 + no" ~ 1] Q@] on 


1 
+ [[g(n, + 1)- 1]? -n.7] @=0 = (77) 


These equations are more difficult to work with and when in a 
system, some of the univalent ions have a value of g not equal to 
one and some have a value of g equal to 1; it is easier to let M and 
A be used only for g equal to 1 and let the other ions be pooled 
into the non-permeating species. It is also easier to let ions with 
other valences be treated as non-permeating. Thus, mathematically 
the ionic species which compose the value of X do not necessarily 
have to be non-permeating. Mathematically, X can refer to all the 
individual species in the system excluding the cation used in 
equation (66) and/or the anion used in equation (67). 


Effect of Changing the Net Valence of Non-permeating Ion 


The parameters 7 and g are influenced by changing the valence 
of the non-permeating ion and can be obtained by differentiation of 
these parameters with respect to the valence. For simplicity, it is 
assumed that the valence of the permeating ions have an absolute 
value of one and that ¢,, = q, =1. Also unless otherwise stated, it 
is assumed that the external environment (side 2) remains con- 
stant, and that finally a new steady state results. 

Differentiation of equation (56) with respect to n, gives: 


tudt --(3!) <0 (78) 
A,] Vo 


Thus, whenever n, increases, 7 decreases. From equation (25), it 
can be seen that a decrease in r means that M, is decreased and 
A, is increased. If M, =M, = 0, then differentiation of equation 


(39) with respect to n, results in: 


— <0 (79) 
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Since n, and n, are both positive, the derivative of r with re- 
spect to n, must again be negative. If A, = A, = 0, then differ- 
entiation of equation (41) with respect to n, is: 


( ") 
dr o 
eee ie eee ef 


; (1-n,)? Ge 


Again, the derivative of r with respect to n, is negative since n, 
and n, are both negative. 
In a similar manner, equation (81) results from equation (58): 


X 
d Rites 
dn, Vo 


It follows from equation (81) that if n, is negative and increases, 
then g must increase. This occurs when A, = A, = 0, as can be 
shown specifically by differentiating Condition IV with respect to 
iy: 
d 1l-n 
payne a (82) 
dn, (1-n,)? 
It also follows from equation (81), that if n, is positive and in- 
creases, then g must decrease. Similarly, this occurs when M, = 
M, = 0, as again can be shown by differentiating Condition IV with 
respect to n,: 
dg (1 +n.) 
—=- ra <0 (83) 
; (1+ n,) 
Effect of ‘‘n,’’ on water shifts. If n, changes to n,, due to ad- 
dition of X ,,, then: 


sishik Tend. 84 
X n Cy 


1 te “te 


Substitution of equation (84) into equation (24) gives: 


) (85) 
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Differentiation of equation (85) results in: 


( ) dg, 

NM, .-7n = 

df, _ zi dase 2 (86) 
dna gn, - 7) 


Combining equations (86) and (81) and (59) gives equation (87): 


* “2 ae 
af, : (n7jz N14) Taig A, + Vo, 
dn 
a g 
(- Zany 
ve (87) 
« (li+n,.)A,, +(1-n,.)M,, 


(- Z Yavin) 


Ifn,, is greater thann,, andifn,, is greater than n,, then X,, is 
greater than zero; and if n,, is equal to one, then equation (87) 
simplifies to: 


le 


Be sey eentte 217) akan, (88) 


dn,, 
-(Z)ana-ay 


a 


Thus, as 7,, increases, water enters side 1. In tissue cells, this is 
the condition which occurs when sodium ions enter the intracellular 
phase, which occurs when the sodium pump breaks down (Conway and 
Geoghegan, 1955). Due to some active metabolic process keeping 
the sodium ions effectively out of the cell (Dean, 1941), the 
sodium ions are essentially the non-permeating ions in the external 
environment. If there are no non-permeating ions in the external 
environment, then such a system cannot be in osmotic balance, 
since g cannot equal zero or infinity. The sodium pump thus main- 
tains the osmotic balance and by altering its activity can cause 
water shifts across the cell membranes. With the exception of 
potassium ions, it appears that the other cations normally present, 
such as hydrogen ions (Caldwell, 1958), calcium ions (Gilbert and 
Fenn, 1957) (Hodgkin and Keynes, 1957) and magnesium ions (Gil- 
bert, 1960a) are also actively extruded from the cell. However, 
quantitatively, sodium is the most important cation. Since cell 
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membranes are extremely permeable to water (Robinson, 1960), it 
seems difficult to imagine any type of direct ‘‘active transport” of 
water. Rather, it seems more likely that the ‘‘active transport 
mechanism of water’’ is really the ‘‘active transport mechanism of 
sodium.” 

When acids are introduced into the internal environment, the hy- 
drogen ions neutralize some of the changes on the non-permeating 
internal anions (i.e., buffering action of the internal anions). Thus, 
the valence of the non-permeating internal ion increases upon the 
addition of acids and decreases upon the addition of alkalies. If 
the value of n, changes just due to an acid effect, then equation 
(24) becomes: 


Ja 
ie = Sr (89) 
g 
Differentiation of equation (89) and the use of equation (81) leads to: 
Se She eee (90) 
dn,, g9 4n,, A, yo, 


Thus, if n,, is negative and increases (due to addition of acid), g, 
increases and /, increases; which means that water leaves side 1. 
This is the situation one would expect to occur in tissue cells. 
However, if n,, is positive and increases (due to addition of acid), 
g, decreases and f, decreases; which means that water enters side 
1. This is what happens to red blood cells. 

Application of equations to muscle cells. For muscle cells, M 
refers to potassium ions and A refers to chloride ions. According 
to the Conway-Boyle theory, the g value for both potassium and 
chloride ions is equal to 1 (Boyle and Conway, 1941). This ap- 
pears to be a fairly good approximation (Hodgkin and Horowicz, 
1959) in spite of the evidence showing that g for potassium ions is 
slightly greater than one (Coombs, Eccles, and Fatt, 1955), 
(Shanes, 1958). 

The concentration of species which can reversibly undergo oxida- 
tion-reduction reactions in the extracellular phase is very small 
and therefore the oxidation-reduction potential of the extracellular 
phase (E,.) can easily be changed. Consequently, the value of Exe 
is dependent, according to equations (16) and (17), upon the value 
of E and the oxidation-reduction potential of the intracellular 
phase (EF, ,). 
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A. Effect of adding acid to muscle cells. Addition of acid causes 
n, to increase and r to decrease (equation (78)). Also, cations 
which bind with anionic sites should cause an increase in n,. 
Perhaps, part of the decrease in r due to the cation, acetyl choline, 
might be due to the binding of it to some specific anionic sites 
causing n, to increase. Decrease in r means that M , decreases and 
A, increases. Equation (90) predicts that f, should increase, 
which means that water should leave the cell. This increase in f, 
tends to increase the concentration of substances on side 1. 

According to equation (13) a decrease in r means that for ca- 
tions, C, decreases if g is constant or g increases if C, is con- 
stant. For sodium ions the value of g is less than 1, so any in- 
crease in g means less work is required of the sodium pump. The 
only way for Na, to decrease is for Na ions to leave the cell. This 
effect of sodium exit would tend to decrease n, (equation (88)) and 
cause more water to leave the cell. 

The relationship between the amount of a species added to side 1 
and f, is given by: 


Zz 
Tees See Oe mek eee | (91) 
On RP r E, qd 
a) C,, = amount added to side 1 times the initial 
volume of side 1 


if 4¢ 1 and 2 = 1, then M is equivalent to C and equation (91) 


a 
M, a ae 


Similarly, if feos 1 and z=-1, then A is equivalent to C and 


oaie Ale are (93) 
A, la fy 


For sodium ions which have been considered to be a component of 
X, the following must also be valid from equation (91): 


Na,, ef ‘e| 7) ee i | (94) 
Na, ey fy INa 


becomes: 


equation (91) becomes: 
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Combining equation (94) and equation (24) and letting X,, =Na,, 


results in: 
o-oo) 
f ms g INa (95) 


: X,-Na, 


If it is assumed that g for sodium remains constant, then dif- 
ferentiation of equation (95) with respect to n,, and substitution of 
equations (78) and (81) leads to: 


d A 
i (49) 2) - Nad Vie 


la 


(96) 


Since n,, is negative, f, increases or water leaves side 1 when 
acid is added. 
The actual relationships between the net transfer of ions be- 
tween the compartments are: 
Hoy wa oo , — Na,, 
=(n,, —7,)X,+(m,, - Na, 
a) H,, = the initial volume of side 1 times the net amount 
of hydrogen ions which entered the cell (side 1). 


(97) 


When acid is added, then H,, is greater than zero. The ratio of r, 
to ris less than one according to equation (78), and f, is greater 
than one according to equation (96). Then from equation (92), it is 
obvious that M,, is negative and if it is assumed that q for sodium 
is constant, then similarly from equation (94), Na,, is negative. 
However, from equation (93), it is not possible to determine if 4,, 
is positive or negative. Actually, the value of A,, can be positive 
or negative depending upon the quantitative relationships. If A,, is 
positive, then A buffers. For example, if itis assumed that Na,, = 0, 
then it is possible to obtain the following: 


dA XiaJ9 
7 SD eee ea Ry. (98) 
neg ra JaV%a g 

Since the ratio of g, to g is greater than one when n,, iS negative 
according to equation (81), then ifn,_ is greater than — 1, the de- 
tivative of A,, with respect to n,, is positive (i.e., the value of 
A,, is positive when n,, is increased). As previously mentioned, 
the value of r decreases which means that A, increases. However, 
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at least part of the increase should be due to the water leaving the 
cell. If the increase in A, due to the water shift is greater than the 
increase due to r; then there will actually be a liberation of A 
(chloride ions) from the cell (i.e., A,, will be negative). Even if 
one excludes a water shift, a slight uptake of chloride ions can 
markedly increase the intracellular chloride concentration since 
this concentration is of such small magnitude. For these reasons, 
the magnitude of A_, will not be as large as the absolute value of 


M.,- When qg for sodium is constant, then it can be seen from equa- 


tions (92) and (94) that: 


el (99) 


Since Na, is much less than M,, the absolute value of Na,, is 
much less than the absolute value of M.,- Thus, it would be ex- 
pected that the uptake of hydrogen ions by the cell would be ac- 
companied by a potassium liberation from the cell, some sodium 
liberation, and hardly any change in any direction of the chloride 
ions. This is essentially what Fenn, Rogers, and Ohr (1958) found. 
For the red blood cell, equation (79) shows that when acid is added, 
r decreases which means that A, (chloride and bicarbonate ions) in- 
creases. The increase in the chloride intracellular concentration 
must be due to uptake of chloride ions since the water moves into 
the red blood cell tending to decrease A,. The tissue cell buffers 
principally by the potassium shift, whereas the red blood cell buffers 
principally by the chloride shift. These same workers did not find 
any decrease in muscle water, but this predicted change would be 
difficult to detect. These ionic shifts contribute to the buffering of 
blood plasma in vivo (Gilbert, 1960b). 

It is evident that the membrane potential is dependent upon addi- 
tions of acids and alkalies. Such changes might play an important 
role with respect to several types of physiological phenomena. 
Thus, the action of insulin increases the membrane potential, which 
secondarily causes potassium to enter the cell (Zierler, 1959). 

B. Application of Conway-Boyle Theory. As pointed out before, 
the non-permeating ion in the extracellular phase is sodium wel 
has a valence of 1, son,=1. The valence of the non-permeating 
ion in the intracellular phase approximates - 1 (Boyle and Conway, 
1941), son, = -—1. This corresponds to Condition III and so g % al 
From equation (49), it is seen that h=1. If r= 40 corresponding 
to a membrane potential of about 93 mv., then from equation (50), 


344 DANIEL L. GILBERT 


the value of a=6 = 8%. Actually n, may not equal — 1, but if it 
were close, g might approximate 1 so closely that experiments 
would not be able to detect any deviation from 1. Thus, the mem- 
brane potential is due to the large intracellular anions which are 
non-permeating due to their large size; and due to the extracellular 
sodium ions which are effectively non-permeating due to some 
metabolic process. Both of these factors contribute about equally 
to the membrane potential. 

If potassium chloride is added to the external environment, this 
is the same as changing A, while keeping n, and X, constant. 


Using equation (56), it can be shown that: 
ar X 
ve ae ga rs Geo ed (100) 


dr 
If n, = 1, and n, <0, then —— < 0 and thus an increase in A, de- 


creases r. Since X,, = 0, equation (89) is valid, and by the use of 
equation (58), it is possible to show: 


df, 1 dg, 
qdA, gdA 
aa 9 Waa (101) 
aie ge SA ae [l(in, - 1) + 9,17 - 9,277] 
Ay. oVo, . * es 
d 
The value of = 0 only when g, = 1 and then no water shifts. 
2a 
If n, = 1, then: 


em Soe oe ST oe fa) 102 
dA,, A,,)29Vo, ; rie. 
Thus, when |n,| <1, water leaves the cell and when (ny loa 
water enters the cell. Experimental evidence shows that no water 
transfer exists (Boyle and Conway, 1941), and therefore the value 
of g must be very close to one. 

If potassium is substituted for sodium in the external environ- 
ment, this is the same as keeping S constant, NM, constant, and 
changing X,. Again by using equation (56), one obtains: 

dr 1 X 

; 1 2 

ae a ra Hy — 9 5 Jing adr mglng - 2) (103) 
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Ifn, = —ng, then g = 1, and: 


dr nS 104 
ax, “2A3 ae 
Decreasing X, results in a decreased r as is evident from equation 
(104). The only time r can vary when the parameters Ny, Ng, and g 
are kept constant is when g = 1. When X,, = 0, then: 


4 = ——*_ = —_* 105 
g AgtAcs aoe ( ) 
or 
Xx g 
fp= st 106) 
1 X, g ( 


From this equation, it is possible to show that: 


1. et es é Le 
nee pie Vo, [- alate 1) oe 
Gre ig je nel | (107) 


If g =1, then |n,| = |n,|. Since the valence is constant, g, = 1, 
and equation (107) simplifies to: 
d 1 1 
aetna /h elon VD (108) 
BX ote Xo, xe 


Thus, when X,, decreases, f, decreases and water enters the cell. 
If it is not assumed that g =1, but that n, is still equal to one, 
then equation (107) becomes: 


af, aay, XoaJ9 2 
Sehwag Fe: aJa (4 _ L 109 
ikexies Sip es a say 


From equation (109), it can be seen that when |n,| <1, then f, de- 
creases when X ,,, decreases and water enters the cell, even though 
g does not aposssacily equal one. Indeed, it has been observed that 
under this condition, water does enter the cell (Boyle and Conway, 
1941). 

In summary, it seems that the value of g for muscle is very close 
to (if not identical with) one. 


This work was aided by United States Public Health Grant A-2669. 
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EXPLANATION OF SYMBOLS 
n 
1+ — 
n 
a= 2 
B 
Lee 
g . 
. A =concentrations of a permeating univalent anion 
n 
1+ 4 
a 
1-9 
. 6B =aconstant called the diffusion coefficient 
rps okas 
n 
Cz 
ZL 
ia 
g 
. C =concentration (actually activity should be used instead of 
concentration) 
C., = amount of substance added to side 1 times the initial 
volume of side 1 
eet a2 
ie 
Led 


. D = concentration of permeating substances 

. E' = observed electrical potential 

. E, = oxidation-reduction potential = -E, =E,,-E 
. LE =E,. — E,, = membrane potential 


h2 


E., = electrical potential 
» =H, — E,,, = reduction-oxidation potential 
woj = Standard reduction-oxidation potential of species j 
wj = Teduction-oxidation potential of species j 
RT oxid 
_ n 
ge EY red 
2 


xX xX 
: Gani (ny- 1)" 75 temaa-Z-+4 


2 
A, 2 


21. 
22. 
23. 
24. 
25. 
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G =nXxX 
hee 
Asda, 
h, = hg, 
es: 
qu 


H = concentration of hydrogen ions 
. | = net pressure on opposite side forcing water into same side 
. J = current density 


31. 


32. 


33. 
34. 
35. 
36. 
37. 


38. 


39. 


40. 
41. 
42. 


43. 


44. 
45. 
46. 
47. 
48. 


m = number of electrons liberated in the transformation of red 


to oxid 


M = concentration of permeating univalent cations 


pe In r7q 
r?7g—1 
Serine 
Ye aides 
r7qg—1 
. K = permeability coefficient 
RT 
Fo ee 
I 


n = valence of X 


Na = concentration of sodium ions 
oxid = oxidized form of species j 


m 


P=-nX 
gq = antiln —1 
A; 
V4 nS as 
M , 
tm = mM, 
Q = Wy 
r = antiln 
r 
Cy = = 
A Tas 
ee 


red = reduced form of species 7 
R = gas constant 

s = distance 

S = total concentration 
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49. ¢ = thickness of membrane 

50. T = absolute temperature 

51. uw, = standard chemical potential 

52. w = total potential 

53. V = volume of water 

54. W = potential due to other forces besides electrical and chemi- 
cal ones 


55. W,=We-W,= work poteatas 


56. wow =n,? 


57. 2 = concentration of effective non-permeating substance 
58. X = actual concentration of non-permeating substance 
59. X,, = amount of non-permeating substance added to side 1 times 
the initial volume of side 1 
60. X,, = amount of non-permeating substance added to side 2 times 
the initial volume of side 2 
61. y = flux 
62. y, = net flux going across the barrier from compartment 2 to 
compartment 1 
63. y, = flux going from compartment 2 to compartment 1 
64. y, = flux going from compartment 1 to compartment 2 
65. z = valence of substance 
66. Subscripts 
a) lrefers to side 1 
b) 2 refers to side 2 
c) a refers to steady state after some change has been made 
d) A refers to permeating univalent anions 
e) e refers to increase in quantity 
f) M refers to permeating univalent cations 
g) Na refers to sodium ions 
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A purely information-theoretical approach to the problem of self-repli- 
cation of elementary living units implies that pure chance is the determin- 
ing factor in the formation of the first living unit. The probability of such 
& spontaneous formation can be calculated from the minimum amount of 
information which an organism must possess in order to replicate itself. 
An estimation of this amount of information is made here by two different 
methods. First by a ‘‘paper and pencil experiment’’ which indicates the 
minimum amount of information needed on a printed page in order that 
with given tools the page could be reproduced. Second—by an analytical 
consideration of some hypothetical molecular mechanisms. A general 
method for handling such problems is suggested. On the basis of esti- 
mated information contents it is shown that under most favorable condi- 
tions the probability of a spontaneous generation by pure chance during 
the lifetime of the earth is vanishingly small. It is concluded that dy- 
namic factors, which may reduce tremendously the information content, 
must play a role in the genesis of life on earth. 


1. Introduction. It is a very frequently expressed point of view 
that the fundamental ability of an organism to reproduce itself is 
tied up to a sufficiently large amount of information contained in 
it. The information must be sufficient to permit the organism to 
select from the environment the necessary ‘‘parts’’ and then to put 
these parts together in a proper order. Such considerations apply 
particularly to the self reproduction of a ‘‘living molecule,” if we 
accept the existence of such a thing. One might perhaps remark 
that the self-replication of genes, which are now identified with 
great probability with the DNA molecule, would make the above 
cautious qualification unnecessary. It may be, however, very 
strongly questioned as to whether in such cases we deal with a 
replication of a single molecule rather than with a system consist- 
ing of several molecules. This point has been particularly empha- 
sized by A. Oparin (1957). Our own theoretical studies of possible 
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physical mechanisms of replication of molecules (Rashevsky, 
1959a,b) strongly suggest the validity of a ‘‘systemic’’ point of 
view. 

On the other hand it must be kept in mind that an approach to 
self reproduction, which does not involve any information theory, 
is also possible, A quantum mechanical model of self reproduc- 
tion, referred to above (Rashevsky, 1959b) is an example. The 
process of reproduction may be the result of definite dynamic inter- 
actions, which determine completely the configuration of the sys- 
tem and which make its information content equal to zero. This 
important point shall be discussed further in the course of this 
paper. 

The pure information-theoretical approach and the pure dynamical 
approach represent two extremes. As most extremes, they are 
probably both inadequate and the likely solution of the problem 
is perhaps best found in a combination of both. For theoretical 
purposes, to help us understand better the possible factors in- 
volved, it may, however, be useful to investigate the two extreme 
cases separately, before trying to blend them properly. We should, 
however, keep in mind that such a study is on a purely abstract 
level. We shall discuss first the information-theoretical approach. 

According to the frequently expressed point of view the require- 
ment of an amount of information sufficient for self reproduction of 
a molecule imposes a lower limit on its size, This is the usual 
‘‘explanation’? of why ‘‘living’’ molecules are very large. In a 
previous paper (Rashevsky, 1955) we pointed out that the problem 
of the minimum size of an organism is analogous to the following 
one: given a written language, how large, in terms of the total 
number of characters, must a book printed in that language be in 
order to contain complete information necessary to manufacture 
that book? 

As we remarked at that time, the answer depends on the nature 
of the environment, on what is given for the manufacturing. The 
same situation applies to an organism, An organism which has at 
its disposal already ‘‘prebuilt’’ organic material needs less infor- 
mation than an organism which must multiply by building up its 
own replica from elementary inorganic material, We should there- 
fore limit our problem to that of a minimal size under relatively 


favorable environmental conditions, in which the number of choices 
is small. 
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2. A Paper and Pencil Experiment. It may be worth while to 
seek an answer to our problem by actually trying to exhibit a set 
of written instructions, sufficient for the reproduction with simple 
means. 

Suppose that the following implements are given: 

1) A set of boxes, one carrying a label ‘‘type,’? others not 
labeled. The labeled box contains type characters mixed with 
other pieces of metal. The non-labeled boxes may contain any- 
thing except type. 

2) A holder, the length of which is equal to the width of the 
page on which instructions are printed. The holder consists of a 
slit into which the different type characters can be inserted in 
arbitrary order. To simplify matters (and thus reduce the necessary 
information) we assume that the slit of the holder and the type 
characters are so designed that while the type can be inserted in 
an arbitrary order, all characters will and can be inserted only 
‘fright side up.’? The holder also has a mark indicating the ‘‘upper’? 
side, 

3) An ink pad of sufficient size. 

4) Several sheets of paper. 

Let us now see how long an instruction, printed on a sheet of 
paper, should be, in order that any person in possession of the 
above described implements should be able to reproduce it, after 
reading the instructions but not knowing before that what he is 
supposed to do. Here is a possible example. 

‘“‘From the labeled box pick out a type character which repre- 
sents the first letter of this text. Insert it into the holder at the 
end marked ‘‘left,’* Then look for another character in a labeled 
box until you find one that corresponds to the second letter of 
this text. Insert it into the holder next to the first one. Repeat 
the procedure, picking out now the character which corresponds to 
the third letter of this text and insert it into the holder next to the 
second character. Continue the above process. When an empty 
space occurs in the text, use a corresponding empty character. 
Continue the procedure until all the letters of the first line of this 
text are in the holder. Hold holder with type down and press it in 
this position against the ink pad. Then take a sheet of paper and 
press the inked holder against it, holding the holder parallel to 
the top edge of the sheet of paper at a distance of 4 cm from the 
edge and with the ‘‘up’’ side of the holder turned upwards. Lift 
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holder from paper, dismantle the type and return it to the labeled 
box. 

“Repeat now the whole process described in the above para- 
graph but refer now to the second line of this instruction sheet. 
When all characters contained in the second line of the text are in 
the holder, press again type to ink pad in the same manner as be- 
fore, and then to the paper, holding it, however, 1 cm below the 
first line which appeared on the paper as a result of the first oper- 
ation. Lift holder, dismantle the type, return to labeled box and 
repeat process with the third line of text, printing the latter 1 cm 
below the second line. Repeat procedure until the whole text is 
reproduced,’’ 

Any literate person in possession of the above described im- 
plements and of the instruction sheet, can reproduce it, even though 
he does not know anything about printing of any kind. The amount 
of information contained in the instructions is sufficient for the re- 
production of the instructions under the assumed conditions. The 
total number of characters in our instructions is 1756. Consider- 
ing the redundancy of English as about 50%, the average informa- 
tion content per character is about 2% bits. Thus the total informa- 
tion is about 4300 bits. This is considerably less than the amount 
of information attributed to a DNA molecule. 

If instead of the above-mentioned implements we would not have 
either holder or ink pad but only the materials from which to make 
them, then of course our instructions would have to be much longer, 
indicating first how to make a holder and an ink pad. Most likely 
we would need several pages of instructions in this case. The 
less ‘‘complete’’ the environment, the more versatile the processes 
involved in ‘‘reproduction,’’ the longer will have to be the instruc- 
tions. The above example, however, probably does represent the 
reasonable minimum of information required in reproduction. 

3. Suggestion for a General Method, If we assume, as is fre- 
quently done, that a living molecule represents a word or sentence 
in an alphabet consisting of individual molecular units arranged in 
a definite order, then the following conclusion appears inevitable. 
The information content of this ‘‘molecular text’’ or ‘‘molecular 
language’’ is determined by the number of letters in its alphabet, 
the possible statistical dependence between some of the letters, 
which makes not all sequences equiprobable and results in a re- 
dundancy, and finally on the total length of the message. The 
same factors determine the information content of a message in 
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any human language. Depending on the differences in the first 
two factors, the same amount of information will require messages 
of different length in different languages. But the amount of in- 
formation needed to perform a certain process is characteristic of 
the process and does not depend on the language in which we ex- 
press it. If in order to describe a process in English we need a 
sentence containing 125 bits of information, then no matter in what 
other language, including the ‘‘molecular language,’ we describe 
the process, the sentence describing it must contain 125 bits. 
Hence if we either actually know the physico-chemical process re- 
quired for self-replication, or if we hypothetize a model of such a 
process, then the verbal description of the process, in any human 
language, will give us the amount of information needed to perform 
this process. A complete verbal description of what a living mole- 
cule must do in order to replicate itself will thus give us the 
amount of information which the molecule must contain. Such a 
complete verbal description is impossible at the present stage of 
our knowledge. However, the procedure may be useful even now in 
estimating the lower limits of needed information content, 

4, An Analytical Approach, Let us now make an approximate 
analytical approach to the problem, Consider a molecular chain of 
any length, which consists of individual molecules of m different 
kinds. We thus have an alphabet of m letters. Assume first that 
there are no dynamically preferred combinations of any letters. In 
other words, assume that the energy € of combination of any two 
‘letters’? A, and A, into A, A, is the same for all ¢ and &, and that 
moreover the addition of a letter A, to a chain A, A, .«.A,;..+A 
requires also the energy €, regardless of the structure of the chain. 
This assumption may perhaps hold approximately but is certainly 
unlikely to hold exactly. This being an abstract study, we shall 
make this assumption at first. With this assumption the informa- 
tion content per letter is 


I, = log, m. (1) 


If the molecular chain contains N constituent ‘‘letters,’’ the 
total information of the chain is 


I =N log, m. (2) 
If this molecular chain, which we shall designate as L-chain is 


to reproduce itself in a manner a living organism does, then it must 
pick up its constituents A,, A,... A,, from the environment, But 
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the constituents of which an organism is built are not as a rule 
contained in the environment in a free state. They form parts of 
other molecular complexes, from which they must be so to say 
‘textracted’’ by the action of the organism, The basic process by 
means of which this is done by organisms on the cellular level, is 
digestion. Even on the level of viruses a process isomorphic to 
digestion takes place, namely the lysis of the host. However, 
even if the constituent signs A, were available in a free state, 
they would be mixed with other molecules, just like the type char- 
acters in the box in our example of section 2. Thus each ‘‘molecu- 
lar letter’? A, will have to be chosen from some n, other possibili- 
ties, and this requires, if all possibilities are equiprobable, in 
other words if no dynamic factor enters into play, an amount of 
information 


I, = log, n, (3) 
for the choice of each ‘‘letter.”’ 


Let the total number of constituents A; in the L-chain be M,, so 
that 


>, =N. (4) 
1 


The total amount of information necessary to pick out the M, 
constituents A; is /;. Hence the total amount of information neces- 
sary to pick out all the constituents to build the L-chain is 


m m 
I, = 1, = Ylog, n,. (5) 
1 1 


After the constituents are chosen, they must be properly assem- 
bled. This process which we shall discuss later on, also requires 
information. Hence the total information content / of the L-chain 
is greater than /., or because of (2) and (5): 


N log, m> » log, n;. (6) 
1 


Consider the average value 7, = M,; log, n, over all indices 7. 
Then we have 


>; log, n, = mlog,n,=m7,. (7) 
1 
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mn», 


N> ; 
log, m 


(8) 


More generally, denoting by q, the amount of information needed 
to ‘‘extract’’ the 7th component A, from some other molecular com- 
plex of which it forms a part, and putting 


m 
Py. (9) 
im1 
we have instead of (8) 
a 
Nie (10) 
log, m 


More generally the information contained in the L-chain will be 
used first to break down a number of molecular complexes which 
contain the different A,’s from amongst the ‘‘debris’’ of those com- 
plexes. In this general case we find, by combining the right sides 
of (8) and (10): 


Nig coed Fix. t9). (11) 
2 m 


For a given environment the quantity 7, in (8), or the quantity 7 
in (10) are constant. Since m is a positive integer, we see that 
the fraction m/log, m has its smallest value, namely 1.89, for 
m=8, The actual minimum of m/log, m is at m = e, and its value 
is e log 2 =1.88. For m=2, or m= 4, the value of m/log 4 m is 
equal to 2. For m=5 the smallest value is 2.15; for m = 8, it is 
2.66. Thus from the point of view of ‘‘optimal design’’ (Rashevsky, 
1960) of the L-chain, the molecular alphabet should consist of 2, 
3, 4, or at most 5 letters, to give the shortest possible chain. If 
the DNA molecule is an L-chain, it has a near optimal design. A 
protein chain is removed from the optimum by a factor of almost 5, 
since 20/log, 20 = 4.63. One might be inclined to speculate that 
this is perhaps the reason why the DNA molecule happens to be 
the basic self-replicating molecular unit, the proteins playing an 
intermediate part. As, however, we shall see presently, this does 
not seem to be likely. 

The excess of N log, m over 7, + q represents the information 
needed to assemble the L-chain from its components 4, after they 
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have been picked up. In evaluating this part of information we 
meet with a peculiar paradox. To arrange the chosen components 
A, in a particular order, assuming statistical independence of 
possible arrangements, requires exactly the same information as 
is contained in this particular arrangement. Hence, we must have, 


because of (2): 


N log, m-m(n,, + 9) 21=N log, m (12) 
or 
m(n_ + Q) 
N ~ ——"——- 2N, (13) 
log, m 


an obviously impossible inequality. 

How does it happen, then, that in our example of a self-reproduc- 
tion of a text this difficulty did not occur? The answer is that the 
process of duplicating a text, as described in section 2, does in- 
volve complex dynamical processes. 

The instruction: ‘‘Insert the type-character, which corresponds 
to the second letter of the text next to the first character’? which 
contains about 250 bits of information actually releases a set of 
conditioned muscular movements. Those movements are not sta- 
tistically independent, otherwise the instruction could not have 
been executed. Similarly the instruction: ‘‘Repeat this procedure’? 
involves a number of complex dynamical processes. But from a 
purely information-theoretical point of view, we exclude all such 
dynamical processes. If we do introduce them, then the difficulty 
vanishes, 

Let us first consider an extreme case. Assume that in an L- 
chain: A, A,...A,... each A; has some side link which results 
in attraction of another A, by a given one. Thereis no attraction 
between A, and A, for i#k. In that case all the free A,’s will 
line up alongside the L-chain, and the chain will thus be dupli- 
cated. There must exist an additional dynamic mechanism to 
secure the separation of the two chains after duplication is com- 
pleted. In this case the information required for duplication proper 
is zero, for there is no choice to make in order to put A; say in 
the ath place, if in the L-chain the ath place is already occupied 
by A,. The process is in this case deterministic. A mechanism of 
that type is found in the replication of DNA, since adenine can 
combine only with thymine, and guanine with cytosine. In the 
building of the complementary chain there is no choice and no 
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more information is required than for example in the formation of 
8 molecules out of an aggregate of H atoms. 

Under the assumption made here to illustrate the extreme case, 
any L-chain would replicate itself if all A,’s are available in the 
free state. The information contained in the L-chain would be 
quite immaterial in this case. However, as mentioned above, and 
as has been discussed previously (Rashevsky, 1956) the ‘‘extrac- 
tion’’ of the components A, from other structures, of which they 
form a part, is an essential property of the living. If the L-chain 
is to exhibit such a property, it may need information stored in it 
in order to produce free A,’s from bound ones. In that case the in- 
formation content of the L-chain is essential for reproduction but 
only indirectly, inasmuch as it is essential for the reproduction to 
obtain free A,’s. 

Mechanisms of reproduction proper, that is of arrangement of 
available free components A, in a proper order which are not of 
such a completely deterministic nature as the one discussed above, 
are in principle conceivable. The information /, required for re- 
production proper will not be zero but may be sufficiently reduced 
below the value /=WN log, m of the information content of the 
chain so that the inequality: 


N logs m-m(i,, + VOSS © (14) 


may be satisfied even though (12) can never be. The information 
], may determine some recursive process. This is the case in the 
example of the reproduction of a text, where we have the instruc- 
tion: ‘‘Repeat this procedure.’’ In that case /, is independent of N. 
Inequality (14) then gives us the minimum length of an L-chain, 
that can reproduce itself in an environment which does not contain 
initially any free components A,. We have 


N>1,+m(a, + (15) 


On the other hand it is possible that /, depends on N, increasing 
with the latter. Then the equality sign in (14) still gives us an 
equation for a minimum N. 

Our conclusion is that while an information-theoretical compo- 
nent may be essential in the theory of self-replicating molecular 
systems, a theory of self-replication based solely on information 
theory without limitation of choice by dynamic factors seems to be im- 
possible. This conclusion seems to be closely related with a similar 
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conclusion made by R. Rosen (1959), in which reproduction is con- 
sidered as a mapping. The relation of our argument to the ‘‘Rosen 
paradox’’ is not yet quite clear but a further study of the connec- 
tion between the two points of view is definitely indicated. 

It may be of some interest to make at least a rough estimate of 
N on the basis of (15), assuming that J, is very small. We then 
have 


NS m(R,, + 9). (16) 


Assuming the design of the chain to be near optimal, that is put- 
ting m=2 and considering that each component A, has to be 
chosen from about a hundred other radicals, we haven, ~ log, 100 ~ 
7. As to 7 let us take the average information necessary to ‘‘ex- 
tract’? a component from some particular molecular aggregate. 
This involves choosing the appropriate aggregate from a large 
number of others. Putting this number again arbitrarily at about 
100 we have for that part of the information again the value 7. To 
that value must be added as a second component of g, the informa- 
tion necessary to loosen the specific bonds which hold the desired 
component in the molecular complex. On that score we can say 
nothing. In line, however, with the discussion in section 3, this 
information can hardly be less than the information contained in 
the italicized sentence which is about 150 bits, provided we assume, 
as we explictly do here, that the whole process is purely informa- 
tion-theoretical and does not involve energy relations. Thus we 
find g~ 157,”,, ~ 7. With m = 2 this gives 


N > 328. (17) 


Since each character of our ‘‘molecular’’ alphabet contains one 
bit of information, the total information of the L-chain is about 330 
bits. This is much less than the amount of information of 20,000 
bits ascribed to a DNA molecule (Lederberg, 1960). 

Let us, however, consider what even that relatively small amount 
of information implies in terms of the probability of an L-chain 
having been spontaneously formed at some time. 

If the information of a configuration of any kind is /, this means 
that that particular configuration is one of 2/ equiprobable con- 
figurations. It may, of course, be one of a much larger number of 
configurations which are not equiprobable. Under the assumption 
of absence of energetic considerations, made here, we have, how- 
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ever, equiprobability. Hence the probability of the configuration is 


p=or! (18) 
With / = 330 we have p = 27 88 ~ 197100, What does such a small 
probability imply? It means that out of every 10 1 configurations 
that will be accidentally formed by the A, components when they 
aggregate into various possible combinations only one in every 
10'°° such combinations will possess the self-reproducing property 
of the L-chain. How much time could it take to produce 10!°° molec- 
ular chains from initially free A, components? That depends on the 
concentration of the A, components and on the conditions and rates 
of the polymerization reactions. Let us make the most generous 
possible assumption, namely that the A, components are present in 
a concentration of 1 mol per liter. This gives about n = 10?! mol- 
ecules A, per cm®. Putting for the average collision cross section 
o* of the A, molecules the value 107 1® cm?, and for the thermal 
velocity »~ 10° cm sec™', the number of collisions of each 4, 
molecule with other A; molecules is about co? vn ~ 107. Out of the 
107 molecules per each cm® about 10'® L-chains of desired length 
of a couple of hundred units will be formed. Making the impossible 
assumption that each collision between two A, molecules leads to 
-a dimerization, or formation of A;A,, and that each collision of a 
dimer leads to higher polymers, we find that in about 1077 of a 
second all initial A; molecules will be dimerized, and a complete 
polymerization into L-chains will take about 10°* sec. Assuming 
further that supply of A, molecules is constantly maintained, we 
see that there will be about 10 L-chains formed in 107* sec, or 
1028 chains per second per cm®. If there should be a sizable proba- 
bility of one L-chain to be formed during the 10° years or 10° sec, 
of the history of the earth, the reacting volume will have to be 
such as to produce 10’ L-chains in 10" sec, or 10% L-chains 
per second, With 1028 L-chains per cm® per sec, this corresponds 
to a reacting volume of 10°! cm’. The volume of the earth is about 
8 x 1074 cm®*! 

Thus, even under the most favorable, actually impossible, con- 
ditions, assumed above, if a self-reproducing L-chain must contain 
even as little as 330 bits of specific information, its chances of 
ever having been formed not only on earth but even in the whole 
universe would be infinitesimal. 

Of course much larger molecular chains are constantly formed 
each having a much higher information content than 330 bits. But 
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what we are inquiring into is the formation of a chain with a spe- 
cific information necessary for self-reproduction. 

Because of the exponential relation between probability and in- 
formation content, the probability increases tremendously with de- 
creasing information. Thus if we need an information content of 
only 200 bits, the probability is p = 2° °° ~ 10°. Under the above 
unlikely conditions we find that one self-duplicating L-chain will 
be produced on the average once in every 10° years in a reacting 
volume of 102! cm8, which is still quite impossible. But with / = 
180 bits, p = 27 18° = 10°54, an L-chain will be produced once in a 
billion years with a reacting volume equal to a layer of 0.01 cm 
covering the whole surface of the earth. A slight further reduction 
of the required information will enable us to obtain about one L- 
chain in a billion years under more realistic reaction conditions 
than assumed here. 

The result of the above discussion is that if life has appeared 
on earth on a purely information-theoretical basis, in other words 
if the first ‘‘living alphabet’’ which carried instructions for self- 
reproduction appeared by pure chance, this first ‘‘living alphabet’? 
must have contained less than 150 bits of information. This is 
much less than is usually ascribed to the information content of 
genes or DNA molecules. 

Large molecules with very specific arrangement of components 
may very well be formed spontaneously if a given configuration is 
energetically strongly preferred to other possible configurations. 
If the energy €, of a given configuration is much smaller than the 
energies €,, &,... €y of the other W possible configurations, 
where W may be very large, then by Boltzmann’s equation, the 
probability of the first configuration may be near one, the remain- 
ing probabilities being very small. Such a configuration may con- 
tain a tremendous number of components but its information con- 
tent will be very small. It is important that in the preceding dis- 
cussion of the probability of formation of an L-chain we did not 
have to refer to the number of components in the L-chain but only 
to its information content. On the basis of our discussion we see 
that a chain with a specific structure and consisting of 10,000 com- 
ponents may be formed spontaneously quite frequently but its in- 
formation content in that case will be small. 

Our estimate of the smallest necessary amount of information in 
an L-chain, 330 bits, is based faute de mieua, on very arbitrary 
assumptions. If we assume the existence of energetic factors in 
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the ‘‘extraction’’ of an A; component from the corresponding mo- 
lecular complex, this may easily reduce the value of g by about 
one half, making it 75 bits, and lowering N to about 160, and the 
total information content of the chain to 160 bits. A ‘*spontane- 
ous generation’? of an L-chain could thus become fairly likely 
during the life span of the earth. 

At first glance it would seem that in the light of the foregoing 
discussion an interpretation of genes or DNA molecules as car- 
riers of large amount of information is inadmissible, since such 
structures could not have been formed. This is, however, not so. 
Once an L-chain does spontaneously appear, it can undergo an 
‘‘evolution’’ in the following manner. The L-chains multiply and 
their number increases exponentially. Assuming for the simplest 
L-chain a generation time of about 1 hour, and a sufficient supply of 
‘*food,’’ within a year there would be 287 ~ 102654 such chains, 
an impossible number. The reproduction will stop long before that 
due to exhaustion of food. But the organic world may well become 
filled with L-chains. Shorter chains, carrying about 20 or 30 bits 
of information will be formed relatively frequently. Their informa- 
tion taken alone, will not be sufficient to result in their self- 
replication. But when added to the information of the L-chain, it 
may result in a useful increase of the total information of the latter. 
With a large number of L-chains available, such additions which 
form new larger L-chains that are richer in information content, 
may be formed relatively frequently. The process of self-duplica- 
tion of the L-chains changes drastically the probability of forma- 
tion of larger chains with specific information. The probability of 
formation of a chain with specific information / is 2~/ and is very 
small if /is large. Let / be just small enough to make the appear- 
ance of the event plausible during the life time of the earth. Let 
I, be the information content of a shorter chain, so that J, <</. 
There will then be plenty of those shorter chains available. The 
probability of a direct formation of a chain with information con- 
tent 1+ /,is 27! x gi — 9° “t!) and, according to the assumption 
made about /, too small. If, however, the chains with information 
I multiply and there is a very large number einen, the probable 
production of the longer chains becomes Nr x 2 1) where Ny is 
the number of the L-chains. This may be a fairly large quantity. 

To sum up, a purely information-theoretical interpretation of 
life seems to be unlikely. Even if life originated with the minimal 
possible information, as discussed above, energetic or other dy- 
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namical factors, must play a role in the process of self-reproduction. 
We thus are led to the study of possible molecular-physical mecha- 
nisms which may be essential for the existence of life. In a pre- 
ceding paper (Rashevsky, 1959b) we gave reasons for believing 
that those mechanisms are to be sought in relational rather than 
metric aspects of molecular structures. In a subsequent paper we 
shall suggest a further development of those ideas which suggest 
other reasons for the large size of a living molecule than its infor- 
mation content. 


The author is indebted to Dr. Robert Rosen for a discussion of 
this paper. 

This work has been aided by United States Public Health Service 
Grant RG-5181. 
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The binding energy of a very long molecular chain, composed of different 
classes of molecules, depends in general on the order of the molecules. 
It is shown that under very general conditions there exists for a given 
brutto chemical composition of a chain, a class of chains which is char- 
acterized by a total binding energy which is equal to the total binding en- 
ergy of any other prescribed chain of different composition within the 
limits of unsharpness of the energy level. This establishes a criterion 
for mapping of a class of configurations of long chain molecules on an- 
other class. To the extent that a mapping constitutes a generalized code, 
those results contribute to the theory of molecular codes. 

Applying to our results the results of a previous paper(1959, Bull. Math. 


Biophysics, 21, 309-326), we arrive at the conclusion that the self- 
replication of a living molecule may be the property not of a particular 
structure but of classes of structures. 


In a previous paper (Rashevsky, 1959) we discussed a possible 
approach to phenomena of molecular self-replication and similar 
ones from the point of view of general energy relations, rather than 
from the point of view of some highly specialized mechanism. We 
came to the conclusion that such considerations in any case per- 
haps provide a clue to the understanding of the fact that the mole- 
cules closely connected with the intimate phenomena of life seem 
to require a minimum size which is rather large. Molecules com- 
posed of only a few atoms do not show phenomena of replication. 
The purpose of this paper is to develop the above considerations 
somewhat further and to suggest some other possible relational ap- 
proaches to those problems. Our discussion may also throw some 
rather different light on the now much discussed problem of ‘‘mo- 


lecular codes.’’ 
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Consider a long molecule, consisting of a chain Ch of N units. 
Those units do not need to be atoms. They may themselves be mo- 
lecular complexes, such as radicals. Let them be, as in the previ- 
ous paper (Rashevsky, 1960) /, units A, in the chain, so that 


> = ¥- (1) 


The binding energy €;, between a unit A; and A, will depend on 
both A, and A,. In general, however, €,, will also depend on the 
nature of the units adjacent to either A; or A,. Thus the binding 
energy €,, between A; and A, in the combination A;A;A,A, will be 
slightly different from the value of €;, in the combination A,A;A,4_. 
Those differences will in general be small, although the differences 
of binding energies between A; and A,(i #4) for different values 
of < and & may be of the same order of magnitude as the energies 
themselves. Denote by € the average binding energy per pair of 
units of the whole chain, so that the total binding energy Ey is 
equal for large N, to 


Ej NE. (2) 


Let us permute the different components A; in the chain. If we 
permute two identical components A,, nothing will change energeti- 
cally. If, however, we permute A, with A,, (¢ 4 &), then, because 
of the above discussed effects, € will change. In general such 
changes are very small. The average € also depends on JN, but 
again only very little. The total number W of permutations which 
do in general produce such a change is 

N! 
We = My (3) 


Hence, if only N and the M,’s are prescribed but not their exact 
order, there are Wy, discrete values of € and hence of E,,, because 
of (2). 

Now consider another chain CA, of N, units B; in which not only 
N, but also the number M,; of the components B, is prescribed. Such 
a chain will possess a total binding energy Ey, =N,€, where €, 
depends again on the particular configuration of the chain. Hence 
there will be Wn, discrete energy levels. 

The values € and €, depend not only on the configurations of the 
chains but also on their length so that Ey, = N€ does not imply pro- 
portionality of Ey and N, 
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5 Let the range of possible values of & for any N be from €& to 
&” > &, while the range of possible values of €, for any N, is from 
&, to €/ > &. If N and N, are sufficiently large, Wy and Wy. will 
be very large and we can consider a distribution function f(® ‘such 
that the number of discrete values of € which fall in the interval & 
€+ d€is given by 


Ww f(®) dé. (4) 


We have 
aa be 
fee ai 8) 25 = 1. (5) 
& 


Similarly, for the chain Ch, we have a distribution function f,(,), 
with 


= 
ik f,(&) dé, i (6) 
4 


Now let us consider the intervals (N€’, NE”) = (EX, EX) as func- 
tions of N. The lower bound of the interval is N€’, the upper is 
Ne’’. For the next value of NV, namely N +1, the lower bound is 
(N + 1) €’, the upper (V+1)€”. If E°S N(E’’- €’) or N2 B°/(E"" - 
€’), then as WN increases all successive intervals (Ey, E‘y) will 
have a non-empty intersection. If €/(€’’— €’) £1, this will hap- 
pen for all N. We shall in the following consider situations in 
which either N is sufficiently large to make N 2 €°/(€”’ — &”) or in 
which €’/(€’’- €) <1. 

For simplicity let us first consider that N, is fixed. We can al- 
ways choose WN in such a manner that the intervals (Né’, Né”’) = 
(Ey, Ex) and (N, €, N,€7’) = (En, » EX) have a non-empty intersec- 
tion. Let the lower and upper bounds of the intersection of (Ey, 
EX) and (EX, E’x,) be E, and E,. Consider a value E such that 


E,<E<E,,. (7) 
Then for sufficiently large values of N and N,, as we have assumed 
them to be, for a class of configurations of the chain Ch, which is 
characterized by a total energy between EF and E + dE there exists 
a class of configurations of the chain Ch which is characterized by 
the same energy between E and E + dE, The number of configura- 
tions of Ch, which belong to this class is according to (4): 


ndi))-HGe 0 
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Similarly the number of elements in the corresponding class of 


configurations of Ch, is 
WN, ( E ) dE : 
N h N ; (9) 


1 1 


From the discussion of the previous paper (Rashevsky, 1959) it 
follows that if the energy E of a configuration of CfA differs from 
the energy of a configuration of the chain Ch, by less than 


AE i 10 
ae (10) 


where A is Planck’s constant and T the average life span of the 
system, then the disintegration of that configuration CA into its N 
components A, may lead to a formation of a corresponding configu- 
ration Ch, from its components B,. 

Hence, to each class C(E, E+ AE) of configurations in CA, 


h 
characterized by an energy between E and E + AE =E + = with E 


subject to (7), there corresponds a class C,(E, E + AE) of config- 
urations in CAh,, such that the disintegration of any element of C 
may result in the synthesis of some element of C,. The number of 
elements or configurations in the corresponding classes are: 


W, rpEt+b/c /E 
ee, Es 
N ds A=) dE (11) 
and 
W E+b/t /[E 
aed — ; (12) 
Nive lM 


We may consider the chain Ch, as consisting of approximately 
N/2 components B,, such that each B, is a pair A,A, which ap- 
pears in Ch, In that case the Ch and CA, chains are identical, ex- 
cept that we consider for each of them a different breakdown into 
components. Now a disintegration of a chain CA into its A, com- 
ponents may result in a formation of two such chains from the B, 
components, though a spontaneous formation may not be possible 
for reasons discussed in the previous paper. 

Similar considerations can be applied to the more complicated 
reactions, expressed by equations (20) to (32) of loc. cit. 
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Theoretically the interesting feature of our results is that we 
may deal with self-replication not of specific molecules but of 
classes of molecules characterized in our case by a given binding 
energy. The chemical composition of all elements of a class is 
the same, but the different elements of the class differ from each 
other in their configuration. 

Following R. Rosen (1960) we may consider in principle, instead 
of the energy, any other conserved quantity, the quantum mechani- 
cal operators which have degenerate eigenvalues. It is difficult to 
suggest off hand what those quantities may be. The interesting 
problem that arises is whether a consideration of several such 
quantities simultaneously will restrict or broaden the classes of 
molecules considered here. 

If something akin to the mechanism discussed here actually 
takes place, we may have to arrive at the conclusion that some, if 
not most of our biologically active molecules are only partially 
specified as to their configurations. Such a conclusion may be 
very far reaching. In view of the complexity of the largest organic 
molecules such a view though sharply different from the accepted 
one, does not contradict any known experimental evidence and may 
therefore merit further scrutiny. 

We have here a peculiar mapping not of a given configuration on 
another given one, but of a class of configurations on another 
class of configurations. To the extent that a mapping represents a 
generalized code, and vice versa (Rosen, 1959a, b), the above con- 
sideration may throw some different light also on the problem of 
‘*molecular codes.”’ 

The ratio of the interval AE = / to the total interval of varia- 
tions of Ey or En, is likely to be very small. Putting T arbitrarily 
at 10° sec, and the width of the range of variation of € (or €,) at 
10-12 erg, we have even for a relatively small value of N, namely 
N ~ 100, AE/(Ex}.- Ey) ~ 10?*. The value of Wy and Wy depends 
according to (3) on the M,’s. It is, however, readily seen that with 
N ~ 100 the value of W, may be very large. Assuming for illustra- 
tion only that we have m different types of components A,, all in 
the same amount, equation (3) becomes: 


N! 


1 eee (13) 
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Using Stirling’s formula, log N = N (log N — 1) we find from (13): 


Wy =m, (14) 


Even with m = 2 this gives W. ~ 10°°. Thus on the average there 
will be 10° discrete energy levels within each interval AEF. It is 
therefore important to emphasize that our considerations do not 
lead to the requirement of a very large number of components in a 
chain. They show, however, that when that number of components 
is sufficiently large, we may have the phenomena discussed. 


The author is indebted to Dr. Robert Rosen for a discussion of 
this paper and for a correction of an error. 

This work has been aided by United States Public Health Serv- 
ice Grant RG-5181. 
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Making some plausible assumptions about the over-all mechanism of 
food catching and consumption by fishes and evaluating in the light of 
those assumptions some available experimental data, it is possible to 
calculate from those data the variation of several important factors with 
the concentration of food. The factors considered are: total rate of 
metabolism, total diurnal energy expenditure in the process of feeding, 
average number of hours per day during which the fish feeds, average 
length of path traveled by a fish per day, and the so-called ‘‘energetic 
coefficient of growth.’’ A possible relation with the work of N. Rashevsky 
(Bull, Math. Biophysics, 20, 299-308, 1959) is discussed. 


The problem of biological productivity of different basins con- 
tinues to be of importance in contemporary hydrobiology. The num- 
ber of investigations devoted to different aspects of this problem is 
rapidly increasing. Those investigations are clearly divided in two 
main groups. The first group comprises the analysis of the processes 
of formation of organic material in the basins (primary production) 
as well as the study of the further fate of the reserves of free 
energy contained in that organic material, as that energy passes 
through a series of heterotrophic organisms. Hence, a realistic 
analysis of energetic transformations reduces to the determination 
of quantitative relations during the transition of energy from one 
level to the next one. In other words we must study the relations 
involved in the flow of energy from the population of organisms 
which in the particular link of the chain of events are used as food 
to the populations of predators. 

It has been shown previously (Ivlev, 1945) that the flow of 
energy from one level to another should be considered from two 
points of view: the ecological and the physiological one. The first 
aspect includes the problems of consumption of food from the sur- 
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rounding medium. The second deals with the analysis of the ener- 
getic transformations of nutrient material inside the organism. 
Other corresponding investigations and generalizations also fol- 
lowed one of the two above-mentioned directions. The ecological 
side of the problem of utilization of food was in particular dis- 
cussed in the books by A. Shorygin (1952) and V. S. Ivlev (1955). 
The physiological side of the problem was discussed in the books 
by G. S. Karzinkin (1952) and G. G. Vinberg (1956). An analysis 
of the basic investigations devoted to a concrete study of the 
production process in different basins is given in a recent revue by 
Vinberg (1959). 

In the present paper the attempt is made to evaluate quantita- 
tively the utilization of food in one particular link of the chain, 
namely fish-plankton. Our problem consisted in the determination 
of the necessary parameters, both of an ecological and of a physio- 
logical nature. We had at our disposal only a very limited factual 
material. This required the use of data which belonged to non- 
identical cases. It also required a number of assumptions which 
are specified below. Thus the aim of the following presentation re- 
duces to the determination of general relations and to the determi- 
nation of the expected orders of magnitudes. We may assume that 
the model of utilization of food presented here can be used, with 
some corrections, in the study of the other links of the production 
process, including the utilization of the primary production. 

As the main object of our analysis we used young specimens 
of Alburnus alburnus (L). A natural population of this species was 
accidentally found in one of the ponds of the hatchery ‘‘Tome’’ 
(Latvian SSR). This population, of unknown size, consisted of 
fishes of the same age and of very uniform size. Other fishes were 
absent from that pond. During the period of observation, from July 
1st to July 18th, at a temperature ¢ = 18.8 - 21.6°C the average 
weight of the young fishes increased from 0.19 gm to 0.81 gm. Thus 
the average weight of the young fish during this period was 0.25 
gm, the total increase of mass of the individual was 0.12 gm, and 
the average daily increase—0.0067 gm. The dry weight of the 
young fishes is about 15% of the total weight. Assuming the 
caloric content of 1 gm of dry material to be 5000 cal (Vinberg, 
1956), we find for the average daily increase expressed in calories 
Q’ = 5.025 cal. 

The young fish in the pond were feeding almost exclusively on 
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Cyclops sp., using both adult specimens and copepodite stages. 
During the period of observation the concentration of the prey re- 
mained approximately constant. The average value was 73455 in- 
dividuals m~°. It varied within the limits from 70040 to 77130 
ind m °. The average live weight of a single Cyclops was 0.015 
mg. Hence the average concentration of the plankton was 1.102 
gm m ® of live weight. Taking into consideration that the ratio of 
dry to live weight of Cyclops is equal to 0.14 (Karzinkin, 1952) 
and that the caloric content of 1 gm of dry material of the fresh 
water copepod plankton is 5140 cal (Ivlev, 1944), we find that the 
concentration p of the food in terms of calories was 808 cal m3. 
Hence the caloric content of one specimen was 0.011 cal. 

It was not necessary to determine the number of feeding young 
fishes because the interaction of fishes of a given population, as 
in the case of any other group of animals of same species and same 
size, can be produced only through the variation of concentration 
of food. Specific intraspecies relations (‘‘complicated’’ competition 
according to Ivlev, 1955) could not have had any important influ- 
ence in the present situation. 

According to our previous investigations (Ivlev, 1945a; 1955) the 
relation between the intensity of nutrition of fishes and the con- 
centration of food is expressed by the following equation: 

dr R 

dp a a ( ci r) o) 
where r is the rate of consumption for a given concentration, FR is 
the maximal rate of consumption for infinite concentrations, and o 
is a coefficient of proportionality. Integrating we obtain 


r=R(I-e %). (1) 


In our study it was necessary to determine empirically the quanti- 
ties R and a. By means of experiments which are described in the 
above-mentioned publications it was found that R = 15.5 calh ‘ and 
a = 0.00037. The value of & was determined by measuring the rate 
of consumption r for a given concentration. It was found that for 
p = 2550 cal m™°, r=9.44 cal h”*. Both coefficients were de- 
termined for a feeding period of 1 hour. 

It is known that the rate of consumption depends not only on the 
average concentration of the food but also on the character of its 
distribution, as expressed quantitatively by the index of aggrega- 
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tion (Ivlev, 1945a; 1955). Inasmuch as we have no knowledge of 
the manner of distribution of the copepods in the pond, we assumed 
it to be uniform. Such an assumption is permissible for copepods 
because as a rule they do not show any tendency to form clusters 
such as found in many species in Cladocera. 

During the process of feeding most fishes, in particular the 
planktophages, move with a definite speed. For the following cal- 
culations we had to determine the speed of locomotion during the 
feeding, and the path traveled during a specified interval of time. 

We have shown previously (Ivlev, 1944), that the velocity of lo- 
comotion of young fishes of a given species and size remains ap- 
preciably constant even when the concentration of food is varied 
within wide limits. The speed decreases only at very high values 
of the concentration p. Using the method described in the previous 
publication we determined the velocity of locomotion of the fishes 
at p = 1000 cal m“° for young Aldurnus alburnus weighing 0.26 gm. 
This velocity v was found at a temperature ¢ = 20°C to be 2.98 cm 
sec’! or 107.28 mh™!. Those and other experiments were used 
with young Aldurnus alburnus taken from the ponds of the Central 
Experimental Station ‘‘Ropsha’’ of VNIORH. 

The above-mentioned parameters refer to the positive part of the 
energy balance. The expenditure of energy which is contained in 
the food is expressed by the values of the metabolic exchange, 
which is calculated from the rates of respiration. 

As is known, the total metabolism is composed of two parts: the 
metabolism in the state of relative rest, in which only spontaneous 
movements occur (standard metabolism Q,.), and the metabolism 
caused by active movements of different intensities (active metabo- 
lism). Therefore we must determine the rates of respiration of the fish 
during relative rest as well as during locomotion with the velocity 
of 2.98 cm sec-'. The method of those measurements and the in- 
strumentation are described in detail elsewhere (Ivlev, 1960). The 
results are shown in Table I. The respirometer which we used per- 
mitted us to determine the rate of respiration at three different 
velocities: 0.65, 2.01, and 6.42 cm sec !. At v = 0.65 cm sec™! 
the fishes did not react to the movement of the water. Therefore 
the rate of respiration was determined only for the two higher values 
of v. As we see from Table I we find the following average values: 
standard metabolism Q , = 0.12 ml of O, per hour; total metabolism 
0.25 ml of O, per hour at » = 2.01 cm sec™* and 1.16 ml of O, per 


hour at v = 6.42 cm sec}. 
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In calculating the rate of respiration at the velocity v = 2.98 cm 
sec”! we used an equation found previously (Ivlev, 1960): 


Q=Q,+mv"; (2) 
where Q denotes the total metabolism, Q, — the basal metabolism, 
while m and n are coefficients. A similar relation has been found 
by Fry (1957). From measurements at two different values of v we 


can determine the values of m and n. They are found to be m = 
0.086, and n = 1.81. Introducing this into equation (2) we find: 


Q = 0.12 + 0.086 y':8!. 


From this equation we find the rate of total metabolism for v = 2.98 
cm sec ' to be 0.38 ml of O, per hour. Hence the ratio of the rates 
of total metabolism to that of the basal is 


Q/Q, = 3.15. 


For further discussions it is convenient to express the metabolic 
rates in calories. According to Ivlev (1939) 1 mg of Q, corresponds 
to 4.77 cal. Hence we have 


Q,=0.57 calh™!; Q=181calh™'. 


From the above basic values we may calculate a number of im- 
portant indices. 

During a period of 24 hours the fish is at times in active motion 
hunting for the prey. At other times it is relatively at rest. It is 
natural to assume that the metabolic intensity during the latter 
period is equal to the standard rate of metabolism while during the 
hunting period, during which the fish moves at a rate of 2.98 cm 
sec ', the rate of metabolism is equal to that of the total. Hence 
if during a 24-hour period the fish hunts for @ hours then the total 
diurnal energy expenditure Q_ is given by 


Q-=a(Q - Q,) + 24Q,, (3) 
or substituting the values of Q and Q, found above 
Q. = 1.24a+4 18.68 cal day '. 


Denoting by Q, the total amount of energy supplied by food (diurnal 
ration), by Q’ — the average diurnal increase of the size of the fish 
expressed in calories we have, according to Vinberg (1956): 


0.8 Q,=Q°+Q-. (4) 
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The coefficient 0.8 characterizes the degree of assimilation of the 
food. From (8) and (4) we have 


0.8 Q, = Q'+ a(Q-Q,) + 24Q,. (5) 
On the other hand we have Q, = ar. Hence 
Vp=ak(I—e “F)% (6) 


From (5) and (6) we have 


Q’ + 24 Q. 


ars 0.8R(1-e° %?)-(Q-Q.) s 


Introducing into (7) the numerical values found above for various 
quantities involved we find: 


a= 9.6 h day~! 


Thus we find that for the given concentration of food young 
Alburnus alburnus feeds 9.6 hours during a 24-hour period. Hence 
it consumes diurnally 38.3 cal, obtained from the food, and loses 
through metabolism 25.58 cal, thus gaining 5.025 cal, while 7.7 cal 
are not assimilated. : 

From the above data it follows that the rate of intake of food per 
day under those conditions equals approximately 0.2 of the caloric 
content of the body of the young fish, since that content is equal 
on the average to 187.5 cal. Therefore the energetic coefficient of 
growth of first order (Ivlev and Ivleva, 1948) is equal to 12.1%, 
while the coefficient of second order is 16.4%. We may remark that 
on the basis of the available data it is impossible to separate the 
rate of consumption per day into the maintaining and the producing 
parts (Ivlev, 1947) and to calculate the third order coefficient of 
both. 

However, the available data permit us to evaluate the coefficient 
of efficiency of food consumption (Ivlev, 1945). In the present case 
this coefficient may be considered as the ratio of the prey organisms 
devoured to those actually met. For p = 808 cal m ° and for the 
caloric content of a single organism g = 0.011 cal, the concentration 
of planktons is 73450 ind m™°. Therefore {73450 = 42 ind are met 
along a path of 1 m. Hence when in the process of hunting the 
young fish travels a path of 107.28 m h~!, it will find along this 
stretch 4506 ind. On the other hand it is known that the rate of 
consumption per hour is 3.99 cal or 363 ind. If we consider that 
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the ratio of the number of the prey individuals consumed to those 
met represents the efficiency of consumption of food, then that 
efficiency equals 8.06%. In other words in the present case each 
fish devours approximately one out of 12 prey organisms. It would 
seem from the above that the presence of the remaining 82% of 
‘“‘non-devoured’’ organisms should not have any effect on the rate 
of consumption. Actually, however, when the concentration of food 
is decreased, the rates of consumption also decrease. For ex- 
ample, if we assume that the feeding of our fishes occurs at the 
concentration which corresponds only to the ‘‘devoured’”’ part of 
the food, in other words at p = 0.806 x 808 = 65.125 cal m ® then 
the rate of consumption per hour will be 15.5 (1- e 9-098" x 65.125) = 
0.4 cal and the rate per day is correspondingly 3.64 cal. The energy 
expenditure during this time is about 32 cal. It is clear that such a 
situation would lead to a rapid starvation and hunger death of the 
fish. 

It is possible that the problem of mutual relation in the system 
predator-prey is very complex and cannot be represented by the 
simple formal scheme presented above. For example, it is not easy 
to evaluate the number of food objects which are ‘‘met.’’ By con- 
sidering this quantity to be proportional to the cubic root of the 
concentration we tacitly assumed that food objects which are even 
slightly off the line of motion of fish cannot be considered as 
“‘met.”’ Actually a planktophage fish notices not only food ob- 
jects which are directly in front of it but also such objects which 
are at some distance from its line of motion. Apparently we must 
consider here two different quantities: on one hand the range within 
which food objects will be noticed and the devouring of which re- 
quires a change of the direction of motion of fish. On the other 
hand the small deviations of the food objects from the ‘‘average’’ 
line of motion, when their catching may be effected by a mere move- 
ment of the body of the fish without essential deviations from the 
line of motion. In the first case the process of catching of the food 
objects will be determined by the length L of the path traveled by 
the fish. In the second case the length of the path remains un- 
changed. Clearly the first quantity will always be greater than the 
second. The difference will depend on the sensitivity of the re- 
ceptors which register the presence of the food objects, as well as 
on the size of the fish. In other words as the fish moves it catches 
such organisms which are located not on a geometric line but within 
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a volume represented approximately by a cylinder the length of 
which is equal to the traveled path per unit time and the radius is 
determined by the deviations of the objects from the path, which 
still permit catching the food without changing the direction of the 
path. 

Making some assumptions we shall now try to pass from the 
particular case to a model of a more general type. We feel that such 
a model should represent a system of functional relations between 
the quantities found in the above particular case and the con- 
centration of the food in the external medium. 

In considering the energy balance we shall take as our unit of 
time 24 hours. Since the time does not enter into equation (1) the 
same type of equation will hold for rates of consumptions per day. 
Expressing the rates of consumption in terms of calories we now 
have 


Q,=R,(1- e7%?), (8) 


where #, is the maximal rate of consumption per day and 4, is dif- 
ferent from a. 

Inasmuch as the real rate of consumption per day Q, = ar there- 
fore from (1) and (8) we find 


R,(1- e %1?) 


Ril — 67%?) » (9) 


In other words the number of hunting hours of the fish is represented 
by the ratio of the rate of consumption per day to the rate of con- 


sumption per hour. 
When &, > a, as is always the case, then a decreases with p. 
The smallest value a_,, of ais obtained for p = ~. It is given by 


a =—., (10) 


On the other hand for p —> 0 the hunting time shall have its larg- 


est value oe. Thus 


To find the value of a ., we expand the exponentials in the nu- 
merator and denominator, preserving only the linear terms. We thus 
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Rio, 
a = 
max R a 
Combining this with (10) we see that 


Oy 


max = — Onn (11) 
Introducing (9) into (3) we find 
= 
Q_ = Set (0 -0,) + 0,. (12) 

By putting p = « we obtain the smallest value of Q_, 

min?- = Fain (Y — Q,) + 24Q,; (13) 
while for p = 0 we obtain the largest value 

maxY~ = Fnax (Q — Q,) + 24Q,- (14) 


Empirically determined and calculated data are available for 1 
point namely for p = 808 cal m *. They are Q, = 38.3 cal, Q_ = 
25.58 cal, Q = 5.025 cal, @=9.6 hours day ‘', and a = 0.00037. 
Data for R, and for 4, are not available. We shalluse equation (8) 
to calculate one of those parameters by more or less arbitrarily but 
plausibly determining the other. 

We probably shall not introduce a very great error by assuming 
the maximum value of the energetic coefficient of the first order to 
be 0.4. Numerous empirical determinations of this coefficient for 
young fishes under normal conditions give close values which vary 
around 0.35. We assume slightly larger values because we con- 
sidered the most favorable conditions for nutritional growth, namely 
the maximal possible concentration of food. In other words, we as- 
sume that for p=, k,=Q’/R,=0.4 Hence from equations (4), 
(8), and (12) we have 


Q’ = 0.8Q,-Q_= 
R (1- e %1P) (15) 
0.8R (1 —e7%1P) - 249 -(Q - nee 
and therefore for p = 0: 


Q’ =0.4R, = 0.8R, -24Q,-(Q-Q,) a (16) 
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Introducing the numerical values for Q, Q,, and R we find R, to 
be equal to 42.8 cal per day. This is approximately 16.5% of the 
size of the fish. With the above found value for the R_ we find from 
equation (8): %, = 0.00278. : 

We further obtain from equations (10), (11), (13), and (14) corre- 
spondingly 


@ nin = 2-15 hours day * 


mi 


@ ax = 20.6 hours day™' 


min’- = 17.09 cal day"! 
max?—- = 39.21 cal day? 


(A) 


(B) 


The above calculated values of the parameters are sufficiently 
plausible, which justifies the assumptions made. 

Now we may represent each one of the above parameters as a 
function of the concentration of the food. They are represented in 
graphical form in Figures 1 and 2. 

The most important parameters are shown in Figure 1. The 
principal elements in our model are the intake and output of energy 
per unit time. We have chosen as such a unit 24 hours. The energy 
balance is thus primarily characterized by the two quantities Q, and 
Q_ as well as by their relation in dependence of the quantity p. As 
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FIGURE 1. 
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seen from Figure 1 with increasing p the rates of consumption per 
day increase according to equation (8) very rapidly and reach prac- 
tically the asymptotic value for p = 2000 cal m *. The curve Q_ = 
f(p), which is obtained from equation (12), is somewhat less steep. 
The point of intersection of the two curves determines that concen- 
tration of the food at which the intake and expenditure of energy 
are equal. It is, however, more correct to determine that point of 
balance from the point of intersection of the curve Q- and 0.8 Q, 
(broken line), because about 20% of the energy taken in is not as- 
similated by the organism of the fish and is excreted in the form 
of non-assimilated products. The point in question can be de- 
termined analytically from equations (8) and (12) by setting 0.89 = 
Q_. It may also be easily determined graphically. We find that the 
concentration of the food corresponding to that point is approxi- 
mately 600 cal m °. Hence below this critical concentration the 
young fishes will starve and eventually die. We may remark that 
for any animal in the process of growth we cannot without harm stop 
the latter for any appreciable length of time (Ivlev, 1955). This 
means that the young fish will perish not at the food concentration 
p’ but at somewhat larger concentration which is needed for the 
smallest possible rate of growth. 

The distances along the ordinates between the curve 0.8 Q, and 
Q. (shaded area) represent the rate of daily growth Q’. According 
to equation (15) Q’ is represented as a function of p by the dotted 
line. We see that this curve also reaches its asymptotic value 
rather rapidly. 

An important parameter is the quantity a, the average number of 
hours per day during which the fish feeds. Apparently it is impos- 
sible to obtain this parameter by direct observation, :at least for 
planktophage fishes. 

We have given above the limiting values of a. It was found 
that as p increases to infinity, a@ does not become less than 2.75 
hours per day. It must be assumed that actually this value will be 
greater because in our analysis we assumed a constant velocity of 
locomotion of the fish during the hunting, whereas actually this 
velocity decreases for very large values of p (Ivlev, 1944). Be- 
cause of this circumstance the curve drawn according to equation 
(9), which shows a rather slow decrease at its right end, will ac- 
tually decrease more rapidly than shown in Figure 1. 

Figure 2 shows some other parameters involved in our model. 
The total length L of the path traveled by the fish per day is equal 
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to av. Assuming for v the value used above, namely » = 107.28m h“! 
we find L = 107.28 a m day !. Hence the smallest value of L, at 
p= is 295 m day ', while the largest value, for p =0, is 2210 
mday *. Taking into consideration the above-mentioned decrease 
of velocity of locomotion at high concentrations the smallest value 
of L should be less than the above given value. 
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Figure. 2 shows also the variation of the energetic coefficient of 
growth of first order, k,, which is nothing else but the ratio Q’/Q,. 
The curve shows that the value of this coefficient, which has many 
times been empirically found to be between 0.33 and 0.35, corre- 
sponds to a concentration of the food equal approximately to 2000 
cal m~? or 2.7 gm m“° of live weight. We notice that this quantity 
actually represents the concentration of food, which is frequently 
found under natural conditions.* 


*It should be pointed out that the above calculations are valid only for 
a particular type of prey, namely the copepods. If the prey is of larger 
size then for the same rate of consumption the fish will require a shorter 
path and will spend less energy, because the distance between the food 
particles is given by 1/h/p. Therefore doubling the size of the prey will 
reduce the distance between the prey objects 8 times for the same value 
of p. Possibly this factor explains the tendency of the animals to choose 


the largest possible food objects. 
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All the curves considered above have one common characteristic. 
The represented quantity varies more rapidly for low values of p 
and less rapidly for larger values. It seems that this peculiarity 
follows from a general biological regularity. The character of the 
trophic interrelations has historically developed in such a manner 
that the concentrations of food actually found in nature are as a 
rule sufficient for nutrition and growth on the higher level. Only 
abnormally low values of p result in an obvious decrease of the in- 
tensity of nutrition and in extreme cases in starvation and dying 
out of the animals. Indeed, cases of clearly inhibited rate of growth 
of animals due to insufficient food supply are not frequently found 
in nature. In most known cases enough food is available to secure 
the nutrition and growth. 

We also notice that all the theoretical curves are characterized 
by the presence of asymptotes for p = ~. Practically those asymp- 
totic values are reached for finite values of concentrations and ap- 
parently those concentrations are within the range of the ones ac- 
tually found in nature. 

Before concluding the present analysis we shall discuss still 
another parameter. We have mentioned above that a planktophage 
fish consumes such prey objects which are not exactly located 
along its line of motion but which are rather located within a cer- 
tain volume. 

We shall make the following assumptions. 

During its locomotion the fish can notice prey individuals which 
are within a certain distance from it. If this distance is sufficiently 
large, then the fish will change the direction of its motion in order 
to catch the prey. This affects the length of the total path traveled. 
When the distance between the prey and the fish is small the fish 
may eat the object without deviating from its path. Let us denote 
this small distance by p and let us assume that at a very low con- 
centration of the prey all the prey objects lying within the corre- 
sponding volume will be swallowed. The volume is equal to 7 p? v. 
Therefore the number of N, of the prey objects that may be poten- 
tially swallowed equals 7p? vp/g. Actually the rate of consump- 
tion per unit time is r=N,q. Hence N, =r/g=R(1-e “?)/q. 
The ratio N,/N, represents the efficiency € of food consumption. 
According to the assumptions made € < 1. When p = 0, then dr/dp = 
a in the numerator. At the same time dr/dp = 7p?v for p= 0. 
Hence for p= 0, 4 =mp7v, and therefore as p approaches 0, € 
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must tend to 1. We thus have 
yee es 


BEE (17) 

Taking the limit of this expression for p — 0, setting it equal to 
1 and using the numerical values obtained above we find p = 4.18 
mm. Hence at near zero concentrations the fish consumes all the 
objects met within a cylinder of a radius 4.13 mm. At higher con- 
centrations some of the objects are not swallowed and the efficiency 
of consumption reduces. We may remark that strictly speaking we 
should speak not of a circular cylinder but of a cylindroid with the 
cross section of 53.3 mm?. 

The curve € = f(p) is shown in Figure 2. Its ordinates were 
calculated from the equation (17). As we see € decreases mono- 
tonically from 1 and tends asymptotically to 0. 

We mentioned above that it is useful to discriminate between the 
“‘noticed’’? and the ‘‘met’’ prey individuals. The first category 
comprises all the prey organisms which lie within the range of the 
corresponding receptors of the predator. This range and therefore 
the distance at which the predator notices the prey cannot be de- 
termined without special experiments. The prey noticed by the 
predator may be caught. For this purpose the fish will have to 
travel a certain distance. The prey may, however, sometimes elude 
the predator, inasmuch as some prey possesses a certain degree of 
motility. 

To the second category belong those prey animals which lie 
within the above discussed volume. Potentially each of those prey 
individuals can be caught. This happens when p is near zero. At 
higher concentrations part of the prey is left alone and only a frac- 
tion characterized by the coefficient of nutrition is consumed. 

We shall now illustrate the above discussed determinations of 
different parameters on a different example, namely the feeding of 
young Clupea harengus membras L. in the Gulf of Finland.* A 
number of desired data are missing and we must therefore make 
more or less plausible assumptions. 

The material was collected during the second part of August. 
The average weight of the young fish was 66.5 mg and the average 


*The data on Clupea were supplied by E. Bityukov, the data on the 
plankton by V. N. Grese. We wish to thank sincerely both of them. 
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daily rate of growth 5.1 mg. Assuming that 1 gm of live weight 
corresponds to 1000 cal we find Q = 5.1 cal day ‘3 

As prey basically small specimens of Eurytemora were used. 
During the period of observation the average concentration of food 
was 6550 ind m~° or 455 m gm? of live weight. In thermal units 
the concentration was thus 327 cal m~°. Assuming that the ratio 
of the weight of the young fish to the weight of the food consumed 
per hour is the same as for Alburnus alburnus we find: R = 4.123 
cal day~! and r= 106 cal hour~! for p = 827 cal m™*. Introducing 
those values into equation (1) we find « = 0.0004. 

To calculate the energy expenditure we shall use coefficients 
which are common to all fishes (Vinberg, 1956). Using the equation 


Q, = 0.3 x 0.0665°°* = 0.0343 ml of 0, per hour 


we find the values of @,. In thermal units Q, = 0.164 cal h7* or 
3.94 cal day™', using 4.77 as the value of the oxycaloric coeffi- 
cient. In line with numerous other examples (Vinberg, 1956; Ivlev, 
1959) we assume that the total rate of metabolism is twice the 
standard one, in other words Q = 20_. Hence Q- = 7.88 cal day '. 
Using these values in equation (4) we find: 


Q = 16.2 cal day™. 


From @ = ar = 1.06a we find a = 15 hour day~!. 


During the period of observation (second half of August) the dura- 
tion of daylight at the latitude of the Gulf of Finland is almost ex- 
actly 15 hours. Hence if we considered our assumptions as correct 
we may conclude that young Clupea feeds all the time during the 
day. 

Clearly the above suggested scheme of tropho-energetic proc- 
esses is not the only one. This is shown by a recently pub- 
lished very interesting paper by N. Rashevsky (1959) in which he 
discusses the mathematical aspects of the process of nutrition in 
fishes. To us that paper is of particular importance because 
Rashevsky analyses a number of concepts which we have formu- 
lated previously (Ivlev, 1955) and which we have frequently used 
in the present paper. 

Without going into details of Rashevsky’s theoretical studies we 
shall briefly discuss some of them. 

The basic principle of Rashevsky’s analysis consists in the use 
of fundamental ‘‘biophysical’’ concepts, in other words in an at- 
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tempt to understand concretely the nature of the different parameters. 
We determine those parameters formally, starting with empirical 
data on the relation between rates of consumption of food and the 
concentration of the food. In particular Rashevsky considers the 
rate of consumption of food as a function of time and introduces into 
his analysis two quantities which he considers as basic. They 
are—the velocity of locomotion of the fish and the distance at 
which the fish notices the prey. 

One should completely agree with the above approach of Rashev- 
sky. Above we also found necessary to use those two basic pa- 
rameters. It was possible to determine the velocity of locomotion 
experimentally. The distance, however, at which the fish and the 
prey interact could be calculated only on the basis of certain as- 
sumptions. Rashevsky quite correctly remarks that a ‘‘biophysi- 
cal’’ analysis of this process requires some special experimental 
work. First of all we should determine experimentally the dis- 
tances at which the prey is noticed by the fish. 

Furthermore Rashevsky considers that the feeding of fishes may 
by analyzed as a process of either a non-stationary or of a sta- 
tionary nature. Our analysis refers to the first type. For this 
case Rashevsky considers the exponential equation (1) which we 
use here as correct. He finds an interpretation of the basic pa- 
rameter and finds 


novgt 


R 


(18) 
where o denotes the distance at which the fish notices its prey, 
t denotes the time, and the other symbols have the same meaning 
as above. 

Rashevsky considers that actually for a given concentration p of 
food the fish eats only a fraction of the available food. As we see 
his picture almost completely coincides with our arguments stated 
above. The difference consists in a certain lack of definiteness 
of the concept of “noticed”? prey objects. Inasmuch as Rashevsky 
considers that in the limiting conditions, ‘‘when the stomach of 
the fish is empty’? each noticed food object is consumed, we may 
consider that those objects belong according to our definition to 
the category of ‘‘met’’ objects. We found for the distance of those 
objects from the line of motion of the fish approximately 4.18 mm, 
Introducing our numerical values into equation (18) derived by 
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Rashevsky we find that this distance is approximately 2.6 cm. At 
the present time it is difficult to decide which of the two quanti- 
ties is more correct. It is, however, clear that this parameter is 
one of the basic ones in the analysis of the process of nutrition of 
fishes under natural conditions. 

From Rashevsky’s theories it also follows that the coefficient 
is not a constant quantity but depends on the length of the feeding 
time. Our calculations confirm this result. As we have seen for a 
feeding time of 1 hour o = 0.000368, while for a feeding time of 1 
day « , = 0.00278. The difference is quite pronounced. 

For stationary conditions Rashevsky considers that an ex- 
ponential relation cannot hold and derives a different one: 

(19) 


r 


% 1+ ap 


In this case 7 as a function of p also tends asymptotically to 
R but the asymptotic value is approached much less rapidly, and 
the curve of equation (18) is always below the corresponding curve 
of situation 1. Indeed attempts to use equation (18) for our data 
gave much smaller values for the ratio of consumption as com- 
pared with those calculated from equations (1) and (8). 

We doubt whether ‘‘stationary conditions’’ actually occur in 
process of nutrition of fishes. The conditions of nutrition are not 
stationary because during 24 hours there will be periods of active 
feeding and periods when the fish does not eat anything. By far in 
most cases the nutrition of fishes shows a definite diurnal rhythm. 
Therefore, with rare exceptions, the nutrition of fishes during a 
sufficiently long period of time should be considered as cyclic but 
not as a stationary process. 

We mentioned above that sufficiently reliable calculations could 
be made only if considerably more experimental material were 
available. Therefore we consider the main contribution of this 
paper as follows: 

First of all we obtained within realistic ranges different pa- 
rameters which characterize the utilization of food by planktophage 
fishes and we found different relations between those parameters. 

Secondly, this analysis has clarified the basic problems the so- 
lution of which requires special experiments. 

We hope that the theoretical models of Rashevsky and the one 
used in this paper shall be useful in the development of a program of 
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further investigations in this sufficiently important branch of 
biology. 
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MATHEMATICAL STUDIES OF THE INTERACTION OF 
RESPIRATORY GASES WITH WHOLE BLOOD 
I. 0, ABSORPTION 


S. R. BERNARD* 
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Some progress has been made on the problem of the interaction of 
respiratory gases with whole blood. A practical working model for oxygen 
absorption in and interaction with whole blood is developed by assuming 
that oxygen molecules compete with protons for binding sites on the 
hemoglobin molecule and by invoking the Wyman-Allen (Jour. of Polymer 
Science, 5, 499-518, 1951) hypothesis that two oxygen molecules go on 
the hemoglobin at one time. Extensive tests of this model against satu- 
ration measurements on blood from humans, horses, oxen and sheep are 
made. Values for the equilibrium constants are calculated and compared. 

In addition a second working model has been developed in an attempt 
to explain why Og saturation measurements when expressed as (100 per- 
cent — percent saturation) are an exponential function of oxygen partial 
pressure. Considerations ares make plausible the following expression 


for saturation, [1— 2 eo Vx/h?, */(1 + (1/20) (8’ nhs + n/p’ ))] are presented. 
Here z denotes oxygen tension, 4 denotes hydrogen ion concentration and 
B’ and y are parameters. 


1. Introduction. This is the first of two papers dealing with the 
problem of the interaction of respiratory gases with whole blood. 
The problem can be illustrated in Figures 1 and 2 which present 
the O, saturation curve and the CO, absorption curve, respec- 
tively, for human blood. It can be seen that CO, suppresses the 
absorption of O, (Figure 1) and O, or oxyhemoglobin suppresses 
the absorption of CO, (Figure 2). It is the purpose in these two 
papers to develop a working mathematical model to predict the in- 
teraction of respiratory gases with whole blood and to test the pre- 
dictions against the available experimental data. The first paper 
considers only the problem of developing and testing the model for 

*USPHS Predoctoral Fellow of the National Cancer Institute. 

**Present address: Health Physics Division, Oak Ridge National 
Laboratory, Oak Ridge, Tenn. 
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FIGURE 1. 


oxygen absorption by blood. The second paper considers the prob- 
lem of developing and testing the model for CO, absorption by 
blood and consolidation of the two models for the interaction of 
respiratory gases with whole blood. 

For a number of years the well-known Adair’s (1925) intermediate 
compound hypothesis based on the assumption that four molecules 
of O, combine with a molecule of hemoglobin has been considered 
to be the basic model for O, interaction with blood. Briefly, Adair 
applies the mass action principle to the following equilibria, 


HbO,,_, + O, = HbO,,, isyee. (1) 


Letting S denote fractional saturation, x denote the O, concentra- 
tion, K ; the equilibrium constant for the jth reaction in (1) and fur- 
thermore letting 


i 
A,K! = IT» 1<j<4 


it is shown that 


s (or) Zn ayceey'], 1<j<4 (2) 
j 
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A, and K are parameters termed an interaction constant and a scale 
factor, respectively. The observation that K, is always greater 
than either K,, K, or K, in order that the value of S calculated from 
(2) agree with the experimentally measured value has been of in- 
terest to investigators for quite some time. This reflects an inter- 
action in binding sites leading to an enhancement of affinity of the 
unbound sites for oxygen molecules. In addition to this effect it 
has been known for some time that oxygenation of hemoglobin liber- 
ates protons. This phenomenon termed the Bohr effect explains the 
suppression in absorption noted in Figures 1 and 2. The Bohr ef- 
fect and the interaction effect have been discussed by many authors, 
e.g., Wyman (1948), Wyman and Allen (1951), Roughton (1954), 
Roughton, Otis and Lyster (1955), Edsall (1957) and Edsall and 
Wyman (1958). 

Adair’s theory cannot be said to be adequate as a working model 
for O, absorption because of the absence of the proton dependency 
factor. On the other hand, Wyman and Allen (1951) have suggested 
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that the Bohr effect and the interaction effect are inseparable and 
that two O, molecules go on the hemoglobin at one time. This sug- 
gestion was supported by structural considerations. In this paper 
we shall show that a quantitative working model for oxygen absorp- 
tion, the interaction effect and the Bohr effect can be obtained by 
combining the Wyman-Allen idea together with the assumption that 
O, molecules and protons compete for the hemoglobin binding sites. 
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FIGURE 3. 


In addition we will develop an alternative model in an attempt to 
explain the observation in Figure 3 that the quantity log (percent 
unsaturation), z.e., the quantity log [100 (1 — S)] is a linear function 
of the oxygen partial pressure. 

In the following we present the derivation and results of the tests 
of these models. 

2. Mathematical Studies of O, Absorption. A. First working 
model. Consider the case where O, molecules and protons compete 
for binding sites on the hemoglobin molecule. Let Z denote a mole- 
cule of hemoglobin, X an QO, molecule and H a proton. Let lower 
case letters, z, x and A denote concentrations. The competitive re- 
actions are 


ZXj;_,+X=ZX,, 1S5j<54 (3) 
H.,2+Hel,Z, l1Sksm, (4) 
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In (3) we have Adair’s hypothesis. Let u, denote the concentration 
of the ZX ,, v, the concentration of the H,Z and u;, the concentra 
tion of the H, ZX, molecules. Let T, and T, denote total concen- 
trations of hemoglobin and oxygen, respectively. The totality con- 
dition for hemoglobin is 


Re a lua) +) oa: (6) 
j k i,k 


while that for oxygen is 


T,=e+)  ju,+)) ity (7) 
j i,k 


where 1<5j <4 and 14m and om denotes double summation. 
i,k 

Let K. denote the equilibrium constants for the reactions in (3), 

L, for the reactions in (4), L ix for the left-hand set of reactions in 

(5) and M,, for the right-hand set of reactions in (5). The mass 

action expressions for these reactions are 


u, = aP, ai (8) 
Ui, = Uy Ni hk =, Q;, (10) 
where 
k 
P, a IT K, Ry, - IT L; 
i=1 i=} 
(11) 


From (6), (8), (9) and (10) we obtain 
Papen Reviews ue vs) | (12) 
k j k 
and from (7), (8) and (10) we find 


Ti =2+ a (: + De Nip a) iP, | (13) 


j k 
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Letting 
Py = (i228 Nar|/(+ Do 2h) (14) 
k 
then from (12), (13) and (14) we find 


T. =e nD iP; P, “\/( +) Ps P=!) (15) 


j 


The range of summation in the above expressions extends from 
1<j<4 and 1<k&<m, assuming, of course, that the hemoglobin 
molecule can bind a total of (4 + m) molecules. 

To obtain the expression for S, the fractional saturation, we de- 
fine 


S =(T, - 2)/(4T,) (16) 
and then from (15) and (16) we obtain 


S= aay sey ( ha ee =) (17) 
j j 


which is the desired expression. It can be simplified for practical 
application to the data by introducing the Wyman-Allen (1951) hy- 
pothesis. Their idea is that two O, molecules go on the hemoglobin 
molecule at one time. In our formulation this can be accomplished 
by letting 


kK, — 0, ae 


18) 
kK, — 0, K, >, ( 
such that 

P, =K, =9, P, = K,K, £0; } (i) 


Thus we see from (19), P, and P, remain finite while P, and P, go 
to zero. In addition we let 


R #0, k=2 an 
* = 0, k= 1,8,4,...,m ( ) 
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and 
#0, j= 2, k=1 
Nip =54=0, f=2, b= 2,8..,m (21) 
= 0, 9 = 4, 1l<k<m. 


and from (20), (21) and (14) we find 
P* = é aoe: Nast) (1+ Rh?) 
k 
PS = (1+N,,A)/(1 + Ry h?) 
= (3+ a Nant) [+ ya”) 
k 


eetlslict Roa)" 


Assume f in (22) is sufficiently large such that A raised to the 
highest power is the dominant term and then we find 


Ps = (Ny /Ra)h™72, PZ = (No1/Ry)/h 


(22) 


(23) 
Py = (Nan/Ra)h™*; Pi = (1/R,)/h?. 
Now let 
Po (Na1/Fq) = Pa’ } (24) 
P.(/R,) = P’’ 


then from (17), (19), (23) and (24) we may write 


‘Sie wal iPr cer! / j + Pe a) (25) 
j 


j 


where now j = 2,4. 
We mention here that (25) can be derived by application of mass 
action principles to the reactions 


ZX gig + 2X = 2Xq,+H 
which are less detailed and less instructive too. Because we will 


have need for the equations leading to (25) we have presented the 


more detailed derivation. 
For convenient application of (25) to data we eliminate h by sub- 


stituting S= 4, e =~,,, solve and obtain 


1 
AY = (PY) 5, (26) 
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and letting 
Ga (PS (Ba (27) 


then from (25), (26) and (27) we find 


5 = (1/4) ae an hy 


1+ 4 (a/a,,)? + (e/z,,)* 


The numerical value for z,, can be read from a plot of the (S,z) 
data and the numerical value for & can be determined from another 
pair of values, say S = 3/4, r=2,,. (Actually, a better but less 
convenient method for the determination of parameter values is de- 
scribed in the Appendix. This method permits one to use all of the 
experimental values for S. It is not as convenient because it in- 
volves the manipulation of matrices.) 

To test this theory the data of Ferry and Green (1929) on buffered 
solutions of crystallized hemoglobin of horse blood can be used. 
From their measurements & was estimated to be 0.8. Inserting this 
value and the appropriate z,, into (28) permits calculation of S as 
a function of x. Figures 4, 5, 6 and 7 present graphs of the theo- 
retically determined values and the measured values of S. Theory 
and experiment do not agree in all cases, e.g., the curve corre- 
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sponding to pH = 7.95. However, we believe for practical purposes 
that the disparity is not significant. 
It is also necessary to test equation (26). Note it can be written 


zs -ipH 
a (it ae RS ae (29) 
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Choosing Pj’ = 25 x 107 20 then (29), as can be noted in Figure 8, 
represents the observed z,, values in the pH range of ca. 6.5 < 
pH < 8.88. Note also the measured scale factors (z,,) at low pH 
values do not fall on the predicted line. This may not be significant 
in view of the paucity of the oxygen saturation measurements at 
low pH (Figure 4). From the more general expression (14) it can be 
shown that the saturation curves become independent of pH at low 
and high pH values, and that the scale factors should level off at 
these extremes. The data are not sufficient to test this prediction. 
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This theory has been tested against additional data. We present 
the results of these tests. 

Figure 9 presents a graph of oxygen saturation measurements on 
human blood at constant pH, the data quoted by Roughton (1954), 
and theoretical curves. Note the following; (a) « = 0.96 rather than 
0.8, (b) theory and experiment do not agree at low O, partial pres- 
sures, (c) Pj’ = 0.8 x 107° (see inset) is quite different from the 
25 x 107?° for horse hemoglobin (Figure 8). 

Figures 10, 11 and 12 present the O, saturation measurements of 
Bock et al. (1924) on human blood x constant CQ, tension and 
theoretical curves. The agreement is excellent. P%’ = 0.47 x 107 2° 
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is little different from the value of 0.3 x 10~?° for human blood at 
constant pH (Figure 9). 

Figures 13 and 14 present the O, saturation measurements of 
Roughton (1936) on buffered solutions of purified ox hemoglobin 
held at constant pH but at different temperatures plotted together 
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with the theoretical curves. The temperature dependence of the 
equilibrium constants can be determined from these data. As can 
be noted from (11) and (24) 


the P3’ and P{’/P3’ are composite equilibrium constants. No at- 
tempt to determine the individual values for the L,, 0, and etc., 
will be made. Instead, letting 


Be = RK and) Pi / PLR KR (30) 
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to denote that these are composite equilibrium constants then from 
(27), (29) and (30) we find 
grP4y? 


Ki Ky = &(2,, 1 
Ke Ke (Kt K3)/(a?) . 
Hence, on inserting the values for the %, z,, and pH into (31) the 
equilibrium constants K{’ Ky and K% kK,’ can be determined. These 
values are presented in Figure 15 which shows them plotted on a 


semilog graph versus temperature. They are conveniently and 
closely represented by an exponential function of temperature. It 
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is to be remembered that the Kj’ Ki and K;’K%{ are composite 
equilibrium constants. We have not considered the possible physi- 
cal significance of this exponential representation. 

Figure 16 presents the very accurately determined O, saturation 
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measurements of Roughton e¢ al. (1955) plotted together with the 
theoretical curve. With a least squares procedure Roughton e¢ al. 
determined a best fitting curve employing Adair’s O, saturation 
equation (2) which would probably agree more ecely with the 
measurements than the curve shown here. However, it is believed 
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that the agreement between experiment and the curve calculated 
with this simple equation is sufficient. 

It is of interest to compare values of the equilibrium constants 
for all of the above data but only after they have been normalized 
to 38°C. We employ the data in Figure 15 to normalize these 
values. To illustrate, for measurements made at 19°C and 38°C 
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there is a factor of ten difference. Hence we divide the values of 
the constants determined from O, saturation measurements made at 
19°C by a factor of ten. 

Table I presents the normalized equilibrium constants for the 
different species on which saturation curves were measured. The 
concentration of the hemoglobin solution employed in the experi- 
ment is also presented. Roughton et al. (1955) states that satura- 
tion curves are a function of the concentration of hemoglobin em- 
ployed in the experiment. We have not made any correction for con- 
centration effect. 


TABLE I 


EQUILIBRIUM CONSTANTS FOR DIFFERENT SPECIES 
NORMALIZED TO 38°C, 


Species Cone. of Hd solution KY Kg Fetes Ky K pA: 
Sheep? 3 mgm/100ce Od 230° 6° 0.33 x rr 
Human Lorca = 5.4 > 6.2 

Human tae 7.3 : 7.9 x 
Ox dilute ? 4.3 » 7.4 3 
Horse® 10.0 Y 17.4 de 


jindicates a correction for euperetuce was made. 
’Dimensions are (mole/1.)/(mm Hg)?. 
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As can be noted in this table the values for sheep hemoglobin 
are smallest while the values for horse hemoglobin are largest. 

This concludes the analysis of the data with this model. Next 
we consider the explanation of the observation that the percent un- 
Saturation is an exponential function of oxygen partial pressure. 

B. Second working model. In an effort to gain some insight into 
the reasons for the exponential representation of percent unsatura- 
tion we consider, in an exploratory manner, the case where we have 
more than one species of hemoglobin. The existence of more than 
one kind of hemoglobin molecule has been suggested by Wyman and 
Allen (1951). In recent genetic studies Itano and Robinson (1960) 
suggest as many as four molecular species of hemoglobin. 

The existence of more than one species is also implied in con- 
sidering the case where the four 0, molecules go on the hemo- 
globin molecule in four distinct steps. We have seen in the pre- 
ceding that only two protons are released on binding four O, mole- 
cules. If we assume that one O, can displace one proton from 
hemoglobin then we are naturally led to inquire which two of the 
four reactions were accompanied by proton release? Viewed in 
terms of a random process we could anticipate 4!/[2!(4 — 2)!] = 6 
distinct ways of removing two protons by addition of four O, mole- 
cules. This implies 6 species. One species releases the protons 
on binding the first and second O, molecules. Denote this species 
by (1,2). Another species releases protons on binding the first and 
third O, molecules, i.e., the (1,3) kind. The remaining four species 
are (1,4), (2,3), (2,4) and (8, 4). 

In deriving the expression for saturation it is not very instructive 
to consider the case where all 6 of the species are present because 
the equation becomes extremely complicated. Instead we investi- 
gate the case where we have the two species (1, 3) and (2, 4) present. 
This will yield some insight into the more complicated case. 

Let Z, denote the ith species of hemoglobin. For purposes of 
completeness and convenience we write the reactions 


Z,X,.,+X=2Z,X, + 41+ (-1C ln (32) 


where 1<j <4 and 1<z<2. When =1, protons are released on 
binding the first and third O, molecules, z.e., the (1,3) kind of 
hemoglobin; and when i = 2, protons are xeleased on binding the 
2nd and 4th O, molecules. Let K;; denote the equilibrium constant 
for the ith species of hemoglobin and the jth reaction in (32), wu ij 
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the concentration of the Z,X; and 
j . . 
ebsites 1 Ses ae See (33) 


Let T. denote total concentration of the ith kind of hemoglobin and 
L 
let 


1 ee ee (34) 
We also have 
4 
Tp = 2; + >a tf (35) 
j=1 
and 
Ti =2+ Mj (36) 
i,f 


for the total hemoglobin and total oxygen concentration, respec- 
tively. The mass action expressions are 


Ui; = 2:9;; (e/htY! (37) 
where 
Cee Aye DIE 1)/) (38) 
From (35), (36) and (37) we obtain 
YA, Qj (0/0 
le a eee (39) 


aes 
i Q; ; (e/h? ¥ 
i 
which as can be easily verified is the same as 


r, a pa is P Q 1j 2a BaP 0) (2/0*Y] z| (40) 
7 j 


Here 0 Sj <4 and Q,,= 1. We now make some plausible statisti- 
cal assumptions in order to evaluate the equilibrium constants. As- 
sume the reactions in (32) in which an X is bound but no proton is 
released, are equivalent and independent, and characterize these 
reactions by a microscopic equilibrium constant denoted by a. 
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Make the same assumptions for the reactions in which protons are 
released but characterize these by a microscopic equilibrium con- 
stant denoted by 6. These constants must be weighted. Assume 
the factors which favor the binding of an O, molecule and the re- 
lease of a proton to be a linear combination of the number (5 — 7) of 
unoccupied QO, binding sites and the average number of bound pro- 
tons which are available for displacement by O,- For the (1,8) 
kind of hemoglobin there is an average of one bound proton and two 
for the (2,4) kind, z.e., the (1,3) kind (first kind) has an average of 
only one bound proton because reduced hemoglobin, Z,, has 2 
bound protons while Z,X,, Z,X, and Z,X, have 1, 1 and 0 bound 
protons, respectively. This makes a total of four bound protons on 
four complexes or an average of 1 per complex available for dis- 
placement. Similar reasoning can be applied to the (2, 4) kind. De- 
note the weighting factor for reactions involving proton displace- 
ment and O, binding by w;;. From the above we set 


w,, = c(5 — 7) + d(2) ; (41) 
where c andd are constants and we assume 
c=d=l1. 


For the reverse reactions, z.e., proton binding and QO, displace- 
ment, denote the weighting factors by w_;, and assume it to be a 
linear combination of the number of bound oxygens, j, and the aver- 
age number of unoccupied proton binding sites, (3 -7), ¢.e., 2 for 


the first kind (1,3) of hemoglobin and 1 for the second kind. Hence 
w_,, = e(j) + (3 - 8) (42) 
where e and f are constants and we assume 
e=f=1. 


These expressions can be used to assign weights to the a and b 
and we thereby determine values for the equilibrium constants for 
the reactions in (32). Since they are equal to the ratios of the 
velocity constants and these in turn are proportional to weights 
where the constant of proportionality is the a or 6 we have 

4 (i= aft) s+ nit!) 
K ,, = (w,,/w_,;,) (a)? (b)? (43) 
With this equation and (33), (38), (41) and (42) specific numerical 
values for the Q,, can be determined. Inserting these values into 
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(40), letting T, =T, performing the indicated multiplication and 
differentiation we obtain 


S_ + (8/3)S” + (87 */8) 8” 
Ti=2+4T on OO (44) 
S, + (B/3)S, + (B ) Sy 
where 
S_= (35/12) a’? + 142°4 + (85/4) a + 2” 
S’ = (5/8) a’ + (42/4) a°% + (75/4) a> + (63/8) 27 
Sy = (9/8)a" + (45/4)a°* + (70/4) 2"* + (85/8) 2" | Oy. 
S, =1+ (85/3) a? + 23274 + (35/8) 2’% + @’® 
So 5a" + 28273 + 80275492" 
Sj =92' + 80079 + 28°54 5277 
and 
pRahes = = oft 
B = (b/a)’ /R?, yy’ = (ad)? » @ = yx/h?, (46) 
Equations (45) can be rearranged to read 
(252°)? (4.22) — (4.32’)® (3.8 27)8 
. +663 0’) (4.8.07 4.52)" “(4.7 9}* 
oi snide io Ga dialiee hie ea 
(1.12) (407) (4.42°)5 (4.327)7 
lg ere oe rears TE, errant Serre 
i} 3! 5! 7! ) 
(47 
: (4.927)? (4.82’)* (4.627)® (3.827)8 
ee ey 20s Bb eee 
(Sa) (5.8w ym (es > tees 
Ss? = oo oO 
i! 3! 5! T! 


sg” (9@’) (5.62’)? (5.12’)5 (4.42°)7 


These expressions are very nearly Taylor series expansions to four 

terms of hyperbolic functions. We write these series in (47) as 
8, = cosh (42’)-1-/f, (2’), 

S, = sinh (42’) - f’(2’), S, = sinh (42) - f/’(2’), 

S, = cosh (4 #’)-f,(2#’), 

S% = sinh (4@’) - fj (@’), SY = sinh (42’) - fy (@’). 


(48) 
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Here the /’s are polynomial residuals, i.e., the difference between 
the given hyperbolic function expanded to four terms and the given 
S polynomial. (Note that the prime symbol on the f does not mean 
differentiation with respect to 2°.) Choosing a value for x’ such 
that the residual polynomials are negligible relative to the hyper- 
bolic functions then from (44) and (48) we obtain 

T,-2 cosh (42’) + (1/8)(8 + 1/8) sinh (42’) - 1 


Expressing the hyperbolic functions in terms of exponentials we 
find 
[1+ (1/3)(8 + B')le**’+ [1 - (1/8)(B + Bole" 47-2 
[1+ (1/8)(B + Bo le 4*’+ [1 - (1/8) (B + Bo) e4*” 


4 


(50) 


» = . oe = + , 
and since e~** becomes and remains neglibible relative to e*4* 
for some value of 2, call it x’’, then 


B= 1—Ye—** 11+ (1/3) (6 + 1/p)l, 2 > 2" (51) 
Letting 
ya4ty,  p’=(b/a)’ ond Un1-S (52) 
and 
A =2[1 + (1/3) (B’/A? + /B NT? (53) 
B = y/BF (54) 


then (51) reads 
dee (55) 


Equation (55) predicts the exponential dependence of U, the 
fractional unsaturation, on oxygen partial pressure. To test this 
model against the data in Figure 3 rearrange (53) and (54) to read 


2/A -1= (B10 FP 4 pr 108? #)/3 (56) 
Goa 0EE (57) 


From a plot of the observed values of B in Figure 3 against the pH 
(which was estimated with the Henderson-Haselbach equation) we 
observe y = 10-5. When we plot the (2/A — 1) employing observed 
values for the A from Figure 3 versus pH we observe a discrepancy 
in agreement between theoretical predictions and observations 
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which is not serious. The right-hand side of (56) is a function hav- 
ing a minimum for some value of A, call ith. When h=h,, then 
Bia= he and the function takes on the value 2/3. This is about a 
factor of 6 or 7 greater than the observed values. In order for this 
theory to fit the data we have to replace the factor of 1/3 by 1/20. 
Although this maneuver would appear to be quite arbitrary it can be 
justified on the basis that the factor of 1/38 is arbitrary, z.e., it 
would be different if we assumed different values for the constants 
c, d, e and f appearing in (41) and (42). From this same plot we 
observe h_ = 107%%. Thus f’ = + 10*?-®. Since we are only in- 
terested in positive values for 4 and B’ we choose B’ = + 10788, 
Both of these values of B’ are mathematically acceptable but only 
one is physically acceptable. Our choice is governed by whether 
we believe a is greater or less than 6 since f’ = (b/a)?- 

To determine the agreement between theoretical predictions made 
with this model and the data we substituted the values for 8’ and 
y into (50), changed the value of 1/3 to 1/20, calculated and 
plotted (50) together with the data shown in Figure 1. Figure 17 
presents this graph. Remarkable agreement is noted for 2 generally 
greater than 20 mm Hg. 

Now we discuss the case where we have more than two species 
of hemoglobin to deal with but only in very general terms. We 
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would expect the polynomial series in (45) to have many more terms 
as we increase the number of species. In (40) where we took the 
product of only two polynomials we obtain an 8 term polynomial. 
For six species we would take the product of 6 four term polynomials 
and obtain a 24 term polynomial. This would correspond to a 
twelve term expansion of the two hyperbolic functions in Taylor 
series or some other function involving an exponential term or terms. 


The author acknowledges the comments, criticisms and sugges- 
tions of Professors H. D. Landahl and N. Rashevsky and of Dr. J. 
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R. I. Macey for his assistance in writing this paper. 


Appendix 


To calculate numerical values of the P, in the equation 


Debye! 
ety Al end (58) 
1+) Pai 


from measurements of the S, call these S,, and the a, call these 2,;, 
multiply (58) through by its denominator, transpose and collect 
terms and obtain 


soak esa 
8,= ). P,G/4 - 8) 2; oe (59) 
j 


N is the number of experimental measurements. To solve (59) for 
the P, let C denote a column vector the elements of which are 


CRS AySay7e2,. Sw) 

where the prime denotes a transposed vector and let 
Eea(Pa eee, Pe 

In addition let A be a matrix the elements of which are 


a,, = (/4 - 8,)«} (60) 


414 S. R. BERNARD 


which matrix has 4 columns and N rows. In matrix form (59) be- 
comes 


C=AP. (61) 
Multiply on the left by A’ to obtain 
A’C =(A’A)P 
and since (A’A) is square its inverse is defined, and thus 
P = (4’A)"*(4’C) (62) 


yields the solution for the P vector of unknowns. For the case 
where A is square, 4 x 4, its inverse is defined and 


PSA AYA Oe Ante. (63) 


The calculations called for in (62) or (63) are laborious but can be 
greatly facilitated with the services of a high speed electronic 


computer. Some more details on this method can be found in House- 
holder (1958). 
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ON THE THEORY OF INDICATOR-DILUTION METHODS 
UNDER VARYING BLOOD-FLOW CONDITIONS 


HERBERT SHERMAN* 
LINCOLN LABORATORY, ** 
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 


The theory of measurement of flow and volume by indicator dilution 
techniques is given in conditions of time-variable flow rates. It is shown 
that the usual Hamilton (1932, Am. J. Physiol., 99, 584-551) methods can 
be misleading if the flow changes at a rate of close to that of the transport 
function. 


I. Introduction. Indicator dilution techniques have been used for 
many years for the measurement of fluid flow (Stewart, 1897; and 
Hamilton, e¢ al., 1932). Theoretical justification has been given for 
this method by Meier and Zierler (1954) and by Stephenson (1948, 
1958). Stephenson, particularly, has used mathematical techniques 
which permit the treatment of the more complex flow problems. 
These techniques will be extended here to the case where flow is 
a function of time. In the human circulatory system, for which 
these techniques are in use, the flow rate varies with the phase of 
the cardiac cycle. While it is intuitively (and empiricaily) evident 
that the flow and volume measured by indicator dilution techniques 
is an average of some kind, it is desirable to have this average 
carefully defined. 

Ii. The transport function. Consider a fictitious capillary bed 
with a single inlet and a single outlet tube. Assume that a quantity 
of indicator M is suddenly added at time ¢ = & at the inlet. Let the 
transport function, f(€, @), describe the rate in which this quantity 
of indicator appears at the outlet; f(€, a) is the fraction of indica- 


*Associate Staff, Peter Bent Brigham Hospital, also. 


**Operated with support from the United States Army, Navy, and Air 
Force. 
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tor which appears at the outlet between ¢ and ¢ + d¢ where 
t=E€+a (1) 


where & = time of indicator addition; a = ¢ — € = age or elapsed time 
since indicator addition. 

Since it will be assumed that all of the indicator added at the 
inlet appears at the outlet* 


mE f(t a, a)ae~ { #(é,¢- é)dt. (2) 


These integrations, which test conservation of indicator, represent 
conditions in which we hold the injection time & fixed and collect 
indicator for all time after injection. 

The transport function may depend on the time** at which the in- 
jection is made with respect to the phase of the cardiac cycle. The 
indicator m,(¢) which appears at the outlet at any time ¢ is a func- 
tion of the entire past history of injection m;(¢) and the variations 
of the transport function f(€, ¢ — €) with time 


mg (t) -{ m,(£) f(é t - 24. (3) 


If the indicator injection is made at time ¢ = €, in a mathemati- 
cally ideal bolus (or impulse function) of the form 


m;(¢) = M8(E - €,) (4) 
where 
um 8(u)du = 1 
5(u)=0; uO 
then 


mo (t) -f Ma(E-£,) f(é,e- dé 
=Mf(é,, ¢ - €,). (5) 


*Even in the case of true laminar flow with zero wall velocity, this 
condition will be satisfied, although the mean transit time may be infinite. 
**The concept of time-varying impulsive response of a system has been 
used elsewhere (Zadeh, 1950) to describe the behavior of linear systems. 
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The resultant curve of the rate of indicator output (not concentra- 
tion) as a function of time is the transport function defined for the 
particular time of injection. 

Consider another indicator injection as a long bolus extending 
from time ¢, to time €, of total quantity Q, and being zero else- 
where as shown in Figure 1. 

The indicator output will then be 


Qh 
m(t) = Teh H(é,¢- 4é. (6) 


This will be recognized as a moving window averaging of f(é, a) 
with a uniform window of length (€, - é,). 


TIME 


INDICATOR BOLUS 


FIGURE 1. 


Any other variation of indicator input as a function of time will 
be recognized as another special weighted averaging of the trans- 
port function. Alternatively and profitably the output may be viewed 
as the distortion of the indicator input variation with time a3 a re- 
sult of weighted averaging by the transport function. 

As a caution it should be noted that the transport function is de- 
fined between two points in a flow system. It may not be, and 
usually is not, representative of the transport function between any 
other two points in the system. 

Ill. The transport function and flow. If the indicator is attached 
to the fluid so that it cannot slip,* then there is a well-defined re- 
lationship between the flow past two points in the system and the 
transport function between these two points. The instantaneous 


*A situation in which this premise is violated occurs in measuring 
plasma flow with indicator-tagged red ceils. The mean velocity of the 
red cells may differ markedly from the mean velocity of the plasma. 
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flow at the output, /,(¢) at a time ¢ is: 
Pod=f FOG e- Oa. (7) 


It is important to recognize that this particular formulation defines 
the output flow in terms of the movement of indicator particles in- 
troduced by F,(€) which particles later appear in F,(¢). With such 
a formulation, one can never define the output flow at time ¢, F, (¢), 
in terms of the input flow at the same time t, unless the channel 
was infinitesimally short or the velocity infinitely high. More 
generally F)(¢) is not measurable by indicator techniques in terms 
of the input flow within an uncertainty equal to the amount of in- 
dicator in transit. 

The transport function f(€, a) and the flow are evidently closely 
related since if the output flow F,(¢) goes to zero then, by the re- 
quirement that there is no indicator slippage, f(&€, a) must also go 
to zero at the same time since there is no indicator output. This 
point will arise again later in a discussion of indicator concentra- 
tion curves. The detailed functional relationship between flow and 
transport function may be determined by an exact knowledge of the 
flow conditions between input and output. Where these flow con- 
ditions are not known, more general side conditions may arise as a 
result of restraints on input and output flow. 

For example, if the flow at input and output is constant as in an 
incompressible fluid flowing at a constant rate through a rigid tube, 
then 


F,(t) = F,(t) =F 


and the transport function must satisfy 


1-[ f(& t- dé. (8) 


IV. Flow from indicator dilution measurements. The usual flow 
measurement is made by adding an impulse of indicator M S(t - a4 
at time ¢= €, and measuring the concentration at the output, the 
output concentration will be: 


ifs Md(é - &,) f(& t- ede 


Mo (t) =oo 


F(t) F, (¢) 


Co (¢) = (9) 
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It is usually assumed (Hamilton, e¢ al., 1932; Meier and Zierler, 
1954; Stephenson, 1958) that there is some constant flow F such that 


(10) 


i] Fc, (dt = M (11) 
and that this constant flow is determinable from 

= M 

| a eee (12) 


fe c, (t)dt 


From the previous equations this calculated flow F is 
— 1 
Ve (13) 
iE f(E4,%- €)) 
Sree ee 
Se Fy (¢) 


where the output concentration must be zero until injection is made 
at ¢= €,. The expression appears to be poorly defined when F, (¢) 
is zero, but if we recall that f(€,, ¢ - €,) must also go to zero at 
this time, then the value of the concentration at such times is the 
concentration in the stagnant fluid which in turn depends on prior 
values of transport functions and flow. 

If the injection is not instantaneous then 


i (14) 
i ad | m,(é)t(& t - dé 
ane lis F, (0 a 


where m,(¢) has been previously defined in Equation (3). 

Several of the limiting cases are relatively easy to study. If 
F,(é) is approximately constant and equal to F,(&,) wherever 
f(é t - €) has a significant value then 


= = = Fy (¢)- (15) 
i) falere i= €,)dt 
gy 
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If F(t) varies periodically with period T which is short compared 
with any epoch in which f(é, ¢ - €) varies, then if Fy (¢) does not 
go through zero, 


~ 1, 1 : 
F = eS esses 
~ fle. C= é,) dt tthe ey) dt 
i F.() ai Fo (0) 
by pup tgee Ota ae 
Ce ee eee 6 
At tt de 
yy et Sileegl ma 
digs teen 
fic tate t-— €,)dt 
&, 
where 
= 1 
nos 1 pe+t de 
al F, (¢) 


In those cases where F(t) goes to zero, these operations must be 
re-examined to find the limiting or stagnant value of f(€,,¢- €,;/F y(t 
and appropriate corrections made in F. 

The rapidly varying periodic flow frequently occurs in the usual 
circulatory measurements since the variations of flow occur at a 
cardiac rate, but the transport function may have significant values 
over several cardiac cycles. 

These equations point out the limitations of flow measurements 
in pulsatile systems. If one deals with flow which varies over a 
time interval having the same order of duration as the transport 
function, then care must be taken in interpreting the flow measured 
by indicator dilution techniques. The experimental indication that 
these variations may be significant is a concentration curve which 
changes appreciably from injection to injection. 

Where the concentration curve varies with the phasing of the in- 
jection, then an average flow, in the sense of Equation (16), can 
be determined by prolonging the input of indicator m,(t) so that 
Fy (¢) changes rapidly with respect to the resultant m,(t). This 
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technique is useful if the prolongation of indicator input does not 
overlap the recirculation of indicator in closed systems like that 
of human circulation. 

V. Volume from indicator dilution measurement. The Hamilton 
method for calculating volume depends on using the concentration 
curve as measured from an impulse of indicator of quantity M added 
at time ¢ = €,, then the Hamilton volume V is defined from* 


V=Ft (17) 


where F has been previously defined in Equation (12) and ¢ is the 
mean transport time. The mean transport time for indicator is 


t(¢,) = tf(€,,t- &,)dt (18) 
1 
(where the integration starts at €, which is the injection time). It 
is evident that this mean transport time is a function of the phasing 
of injection relative to the cardiac cycle. 
Unfortunately the transport function is not usually available 
directly. All that is usually measured is 


oo oe tf(€,, ¢- €,) 
te (t)dt = M —_——_————-qi. 19 
i ae | Fy (2) Ose 


Thus, the mean transport time that is usually measured is ac- 
tually weighted by the variable flow rate and the calculated volume 
is, from (17) and (11) 


ie etetyde 
fase 


V = ——______, (20) 
| f e(oat| 
0 
ie tf(€,,t - €1) dt 

Te LX er 


fie dese t= é)) ae) 
0 F, (2) 


*This definition assumes that the transport function falls faster than 
1/t? in the limit, that a mean transport time exists, that there is no ins 
dicator slippage with respect to fluid, and that the concentration becomes 
everywhere the same in the volume under measurement (Stephenson, 1958). 
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One can make similar remarks, as have been made concerning 
flow, about the mean transport time if the flow changes much more 
rapidly than the transport function or vice versa and the same cau- 
tion should be exercised about the significance of the Hamilton 
volume when the flow changes at rates similar to the transport 
function. 


I am indebted to Mr. William Levison of the Research Laboratory 
of Electronics, Massachusetts Institute of Technology, for his 
perusal and comments on the manuscript and to Dr. Lewis Dexter 
for discussions leading to this problem. 
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RECEIVED 4-12-60 


ANNOUNCEMENT 


SYMPOSIUM ON MATHEMATICAL THEORIES OF 
BIOLOGICAL PHENOMENA 


A Symposium on Mathematical Theories of Biological Phenomena 
will be held at the Barbizon-Plaza Hotel in New York City on May 8, 
9, and 10, 1961 under the joint sponsorship of the Committee on 
Mathematical Biology of the University of Chicago and of the New 
York Academy of Sciences, under the general Chairmanship of Pro- 
fessor N. Rashevsky. 

About 12 one hour lectures will be given on different topics of 
mathematical biology with ample time left for discussions. The 
lectures will be on a technical level. The following is a partial 
list of lecturers and of the subjects treated: 


A. F. Bartholomay, Harvard Medical School—Mathematical Theory 
of Reaction Rates 

L. Danziger and G. L. Elmergreen, Milwaukee Sanitarium Founda- 
tion and University of Wisconsin, Milwaukee—Theory of Peri- 
odic Catatonia 

J. G. Defares, University of Leiden, Holland—Theory of Lung 
Functions 

E. H. Kerner, University of Buffalo—Gibbs Ensembles and Bio- 
logical Ensembles 

H. D. Landahl, University of Chicago—Mathematical Theory of the 
Central Nervous System 

S. Roston, University of Louisville Medical School—Dynamics 
of Cardiovascular Phenomena 

G. Sacher, Argonne National Laboratory, Chicago—Theory of 
Mortality 

J. H. Woodger, University of London—Axiomatization of Biology 
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